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In  this  thesis,  I study  the  conformational  properties  of  polymers  in  ncu-' 
stationary  and  stationary'  states. 

In  the  fiist  part  of  my  thesis,  I study  the  conformations  of  polymers  in  non- 
stationary state.  I consider  two  cases:  (1)  the  motion  of  a polymer  penetrating 
a membrane  and  (2)  the  dynamics  of  DNA  molecules  in  a get.  In  the  first  case, 

I assume  that  a polymer  crosses  a membrane  several  times  during  the  process 
of  penetration.  The  polymer  reptates  along  its  own  contour.  The  mobility  of 
the  polymer  chain  is  related  to  its  configurations  which  are  characterized  by  the 
number  of  crossings  with  the  membrane.  I find  the  configurations  which  give 
the  main  contribution  to  the  transport  and  calculate  the  overall  permeability  in 
the  case  of  both  Zimm  and  Rouse  dynamics.  As  a second  system,  I consider 
the  motion  of  a long  DNA  molecule  confined  in  a gel.  When  a strong  electric 


field  is  applied,  the  electric  field  forces  the  DNA  into  a tree-like  structure.  The 
dynamics  of  large  segments  of  DNA  is  almost  deterministic  and  can  be  described 
by  a set  of  simple  mechanical  equations.  This  allows  the  numerical  study  of  gel 
electrophoresis  of  DNA.  1 confirm  the  hypothesis  of  the  statistical  self-similarity 
of  a moving  polymer. 

In  the  second  part  of  the  thesis,  there  are  three  major  sections  all  dealing 
with  configurational  properties  of  a single  isolated  charged  polymer,  1)  I study 
the  configurational  statistics  of  a ring  polyampholyte  chain.  In  addition  to  the 
cascade  like  transitions  similar  to  that  in  a linear  chain,  I found  that  a ring 
polyampholyte  has  a metastable  phase  where  the  configurations  are  controlled 
by  the  randomness  in  the  charge  distribution.  2).  I study  new  exact  solution 
of  the  nonlinear  Poisson-Boltzmann  equation  describing  counterion  condensation 
on  a strougly  charged  polymer  in  the  dilute  regime.  The  counterion  has  to  be 
brought  from  three-dimensional  space  and  the  entropy  loss  is  3/2  times  larger 
than  in  the  picture  of  Manning.  The  condensation  occurs  at  a critical  line  charge 
density  A = 3/2,  which  is  3/2  times  larger  than  in  Manning’s  theory.  3).  I study 
the  rigidity  of  the  charged  polymer  in  the  presence  of  counterion  condensation, 
I show  that  the  condensation  of  counterions  changes  the  local  structure  of  the 
charged  polymer  more  than  in  the  Debye-Hiickel  approximation.  By  solving  the 
Poisson  Boltzmann  equation  in  cylindrical  coordinates  in  the  vicinity  of  charged 
polymer,  I show  that  the  persistence  length  lp  of  a strongly  charged  polymer  is 
proportional  to  the  Debye  radius  Ip  ~ Rj. 


CHAPTER  1 


INTRODUCTION 

A polymer  is  a macromolecule  consisting  of  many  repetitions  of  the  same 
basic  unit  called  "monomer”  [1-4].  For  example,  polyethylene  is  the  polymer 
made  of  ethylene  molecules  (monomers),  A DNA  molecule  consists  of  up  to  107 
nucleotides.  Polymers  are  very  common  in  our  life.  Plastics,  rubber,  proteins 
and  other  biomolecules  are  polymers.  A material  made  of  long  molecules  has 
very  different  properties  from  one  made  of  small  molecules.  Polymeric  materials 
are  very  flexible  and  viscous,  and  elastic.  However,  because  of  their  complexity, 
polymers  have  been  studied  only  for  a few  decades. 

The  conformation  of  a polymer  in  solution  depends  on  the  solvent,  concen- 
tration and  temperature.  In  a good  solvent,  the  interactions  between  monomers 
and  solvent  are  attractive;  therefore,  polymers  can  dissolve-  The  effective  inter- 
actions between  monomers  are  repulsive  so  that  polymers  stretch  out.  On  the 
other  hand,  in  poor  solvents,  the  net  interactions  between  monomers  is  attrac- 
tive; therefore,  precipitation  occurs.  Depending  on  the  concentration,  a polymer 
solution  can  again  be  classified  into  three  different  regimes:  dilute,  semi-dilute 
and  concentrated  solution. 

Experimentally,  in  order  to  characterize  the  polymer  (for  example,  to  measure 
the  molecular  weight),  one  places  the  isolated  polymer  in  a dilute  solution  where 
the  interactions  with  other  polymer  chains  are  neglected. 


However,  a typical  polymeric  material  is  rubber  like  and  many  polymers  are 
entangled  with  each  other.  In  order  to  understand  the  properties  of  polymeric  ma- 
terials, we  have  to  consider  many  chain  effects.  Polymers  are  long  one  dimensional 
objects  and  they  can  not  cross  each  other.  In  a many  chain  system,  the  topological 
interactions  between  polymers  are  important.  Topological  interactions  strongly 
affect  the  dynamical  properties  because  they  impose  topological  constraints  on 
the  motion  of  polymers.  With  new  concepts  introduced  by  P.  G.  de  Gennes  [1] 
and  S.  F.  Edwards  [2],  the  framework  to  understand  the  many  chain  systems  - 
like  concentrated  solution  and  melts  - has  been  established.  There  are  two  models 
of  the  dynamics  of  unentangled  polymers  - Rouse  and  Zimm  dynamics.  In  the 
Zimm  model,  the  monomers  are  coupled  to  each  other  by  hydrodynamic  inter- 
actions, while  the  Rouse  model  assumes  that  there  is  no  hydrodynamic  coupling 
between  different  monomers  [1,2]. 

In  the  first  part  of  my  thesis,  I study  the  conformational  properties  of  poly- 
mer in  dynamical  processes  in  the  presence  of  topological  obstacles.  In  chapter 
2, 1 study  the  conformational  properties  of  a polymer  penetrating  a membrane. 
There  is  an  interplay  between  dynamic  and  conformational  properties.  The  dy- 
namics of  the  polymers  is  discussed  in  both  the  Rouse  and  the  Zimm  regime. 
In  chapter  3,  the  dynamics  of  DNA  molecules  - ultra  long  polyelectrolytes  - in 
gel  electrophoresis  is  studied.  The  stretched  tree-like  structure  of  the  DNA  con- 
formation under  the  strong-held  electrophoresis  provide  deterministic  equations 
for  their  dynamics.  Through  theoretical  analysis  and  the  numerical  solution  of 
these  coupled  differential  equations,  I calculate  the  average  velocity,  the  period  of 


motion,  tho  loop-size  distribution,  the  critical  velocity  of  the  stability  of  tree-like 
structure  and  other  properties. 

In  the  second  port  of  my  thesis,  I study  the  conformational  properties  of  a 
single  isolated  charged  polymer.  Polymers  consisting  of  monomers  which  include 
ionizable  groups  are  called  “polyelectrolytes."  The  Coulomb  interaction  between 
the  monomers  affects  the  microscopic  structure  of  polymer,  which  results  in  very 
different  properties  of  polyelectrolytes  compared  to  that  of  neutral  polymers. 
The  main  theoretical  difficulties  arise  with  the  long  range  nature  of  the  Coulomb 

This  part  of  the  thesis  (chapters  4 ~ 7)  is  organized  in  the  following  way.  In 
chapter  4,  the  general  difficulties  in  dealing  with  charged  polymers  are  discussed. 
In  chapter  5,  1 discuss  the  conformations  of  a ring  polyampholyte.  The  confor- 
mations of  polyelectrolytes  and  polyampholytes  have  been  studied  by  Dobrynin 
et  al  (5|.  and  Kantor  et  al.  (6,7],  respectively.  A polyampholyte  is  distinguished 
from  a polyelectrolyte  because  it  carries  both  positive  and  negative  ions  along 
its  backbone.  There  are  similarities  between  the  static  configurations  of  a lin- 
ear polyampholyte  and  polyelectrolyte.  After  a short  review  of  these  previous 
works,  I discuss  the  randomness  dominant  phase  found  in  a ring  polyampholyte 
chain.  I use  Monte  Carlo  simulations  in  order  to  illustrate  the  properties  of  the 
conformations.  Simple  mean  field  type  arguments  are  presented. 

The  analytic  approach  to  solve  the  nonlinear  Poisson-Boltzmann  equation  is 
illustrated  in  chapter  6.  A new  solution  of  the  Poisson-Boltzmann  equation  in  the 
dilute  solution  regime  provides  a new  picture  of  the  condensation  of  counterions. 


I compare  the  results  with  the  classical  condensation  theory  [8]  suggested  by 
Manning. 

In  Chapter  7, 1 analytically  study  the  effect  of  the  nonlinearity  of  the  Poisson 
Boltzmann  equation  on  the  persistence  length  of  polyelectrolytes. 

The  results  are  summarized  in  the  conclusions,  and  topics  for  further  study 
are  discussed. 


CHAPTER  2 


DIFFUSION  OF  A POLYMER  THROUGH  A THIN  MEMBRANE 

In  a recent  paper  Yoon  and  Deutsch  [9]  considered  an  interesting  and  im- 
portant problem  of  diffusion  of  a linear  polymer  in  the  presence  of  permeable 
membranes.  This  type  of  transport  is  the  key  element  of  many  biological  and 
technological  processes.  It  has  the  potential  to  be  used  for  the  separation  of 
polymers  of  different  lengths  [9]  and  different  compositions.  It  can  also  be  used 
for  controlled  drug  release.  Usually  a membrane  is  thought  of  as  a thin  materia] 
interface  separating  two  solutions.  Such  a membrane  is  a strongly  fiuctuating 
object  which  might  have  a rich  variety  of  different  phases  by  itself  [10-16].  The 
mechanical  properties  of  the  membrane  might  change  with  the  adsorption  of  a 
polymers  on  its  surface  [17].  Besides  it  is  a mechanically  fragile  object  which  can 
make  it  difficult  to  perform  well  controlled  macroscopic  measurements. 

A simple  way  to  avoid  this  problems  is  to  make  a two  dimensional  “membrane" 
by  putting  together  two  slices  of  gel.  The  surface  of  the  gel  can  be  treated  in  such 
a way.  that  after  two  slices  of  a gel  are  brought  together  they  form  an  internal 
interface.  This  membrane  is  mechanically  very  stable  and  its  permeability  can 
be  controlled  by  various  surface  treatment  techniques. 

I consider  a dilute  solution  of  polymer  chains,  each  made  of  N monomers. 
The  length  of  each  monomer  is  6.  The  bulk  concentration  pa  of  polymers  on 
one  side  of  the  membrane  is  maintained  constant  and  is  zero  on  the  other  side. 
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The  polymers  are  penetrating  the  membrane  which  has  negligible  thickness  and  a 
mesh  size  a of  order  6.  The  membrane  has  a chemical  potential  n measured  with 
respect  to  the  bulk.  The  chemical  potential  will  always  be  given  in  kgT  units.  The 
statistical  weight  of  any  polymer  configuration  with  k crossings  of  the  membrane, 
compared  with  the  similar  configuration  in  the  bulk,  has  an  additional  factor  e~ku 
associated  with  these  crossings.  1 neglect  all  excluded  volume  interactions.  For 
dilute  solutions  1 assume  that  the  bulk  chemical  potential  on  both  sides  of  the 
membrane  is  zero.  Later,  in  order  to  describe  filtration  of  concentrated  polymer 
solutions  through  the  membrane  I assign  different  bulk  chemical  potentials  on 
different  sides  of  the  membrane. 

Yoon  and  Deutsch  [9]  considered  the  simplest  case  when  the  chemical  potential 
on  the  membrane  is  zero  (p  = 0)  and  assumed  that  the  configuration  of  the 
polymer  crossing  the  membrane  is  Gaussian  (see  Fig.  2.1  a).  In  this  “typical" 
configuration  the  polymer  crosses  the  membrane  \/N  times  and  its  mobility  is 
greatly  reduced  compared  to  the  mobility  in  the  solvent.  The  configurations 
with  fewer  crossings  are  thermodynamically  unfavorable,  but  they  have  larger 
mobility.  I will  show  below  that  if  effects  of  hydrodynamic  friction  are  taken 
into  account,  these  configurations  give  the  dominant  contribution  to  the  total 
transport  through  the  membrane.  The  chemical  potential  of  the  membrane  acts 
as  an  additional  control  parameter  in  this  problem.  If  p is  positive  (negative)  it 
will  result  in  an  decrease  (increase)  of  the  number  of  crossings  k and  affect  the 
transport  through  the  membrane. 


Figure  2.1:  Configurations  of  the  polymer  during  penetration  of  a membrane, 
a)  Gaussian  configuration . The  number  of  crossings  k in  this  typical  configuration 
is  k — Af1/S.  b),  c)  Sketch  of  a typical  configuration  when  p > 0 and  p < 0. 
The  number  of  crossings  k is  controlled  by  the  chemical  potential  p.  d)  Effective 
tube  picture.  The  motion  of  a polymer  can  be  thought  of  as  a reptation  through 
a series  of  “gates" . The  distance  between  the  gates  is  of  order  6 (JV/fc)1/2. 


2.1  Reptatiou  Time 


The  mobility  of  a polymer  which  crosses  the  membrane  k times  (see  Fig.  2.1 
d)  can  be  calculated  following  the  analogy  with  rcptation  in  a tube  [9,18,19].  The 
typical  length  of  a polymer  between  the  two  subsequent  crossings  of  a membrane 
is  N/k , and  the  physical  distance  between  the  crossing  points  is  6 jN/k  (see  Fig. 
2.1  d).  The  diffusion  time  through  the  membrane  can  be  calculated  as  the  time 
needed  to  reptate  through  k "gates",  with  a distance  b -J N/k  apart  from  each 
other.  The  total  length  of  this  hypothetical  tube  is  L = b -/Nk  and  the  time 
needed  for  a polymer  to  reptate  through  the  tube  is  i = i!/D  = tfkN/D.  Here, 
D is  the  diffusion  coefficient  of  a polymer. 

The  simplest  approximation  for  D is  given  by  the  Rouse  model  for  mobility 
with  D = DjN  and  reptation  time  1*  = b*kN~ / Dn.  For  a Gaussian  configura- 
tion with  k = -/N  this  will  reproduce  the  Yoon  and  Deutsch  [9]  result  t ~ Nin. 

The  application  of  Rouse  dynamics  can  be  justified  in  special  cases  when 
the  hydrodynamic  interactions  between  monomers  of  the  chain  can  be  neglected. 
This  can  happen  due  to  the  presence  of  other  polymers  (semidilute  solution) . or  if 
the  polymer  is  confined  inside  the  gel  which  damps  the  hydrodynamic  interaction 
between  monomers  of  the  chain.  In  the  latter  case,  the  membrane  is  actually  an 
interface  between  two  pieces  of  a gel  brought  together  as  described  before. 

If  hydrodynamic  effects  are  important,  the  diffusion  coefficient  is  controlled  by 
viscous  friction  at  the  gates  where  the  polymer  penetrates  through  the  membrane 
as  shown  in  Fig.  2.2.  In  order  to  calculate  the  friction  I assume  that  the  polymer 
is  pulled  with  a small  velocity  V through  the  gate.  This  velocity  should  be  small 


Figure  2.2:  Illustration  for  the  calculation  of  hydrodynamic  friction.  The  polymer 
is  pulled  through  the  gate  with  small  velocity  V.  The  friction  on  a polymer 
fragment  of  length  .'V,  is  about  the  same  as  that  of  a sphere  with  radius  R,  — 6.V1  /- 
moving  with  a drag  velocity  of  Vt  = R^V/Ni  ~ 1/R,.  Most  of  the  dissipation  of 
energy  occurs  at  a minimum  scale  ~ a,  (gate  size).  The  contribution  from  large 
scales  is  negligible. 
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enough  to  treat  the  configurations  of  the  polymer  as  if  it  were  in  equilibrium.  The 
friction  on  each  section  of  a polymer  of  length  Nt  (as  counted  from  the  gate)  and 
size  ft,  = -/W,  can  be  calculated  as  the  viscous  friction  of  a sphere  of  radius  ft, 
pulled  with  velocity  Vt  = R,V/N,  towards  the  gate.  The  friction  is  / = 6ir RqV) 
(q  is  solvent  viscosity).  The  rate  of  energy  dissipation  due  to  this  friction  is 
ft  = fiVi  ~ 1/fti.  This  means  that  the  main  contribution  to  friction  comes  from 
the  minimum  scale  ftj  which  is  the  gate  size  o. 

If  there  is  no  large  dispersion  in  the  sizes  of  different  gates,  the  dissipation  of 
energy  at  each  gate  is  about  the  same  and  the  total  dissipation  is  proportional 
to  the  number  of  gates  k.  The  overall  mobility  of  the  polymer  is  ~ 1/jfc  and  the 
diffusion  coefficient  can  be  calculated  from  the  Einstein  relation  between  mobility 
and  diffusion:  D = kBTi/.  The  reptation  time  in  this  case  is  tf  = fik7N/D„. 
Here  Do  is  defined  in  terms  of  solvent  viscosity  q:  D„  = kgT/ri  . I will  see 
that  because  of  the  k‘  dependence  of  the  reptation  time  tf  the  contribution 

suppressed,  and  configurations  with  few  crossings  give  the  main  contribution  to 
the  membrane  transport.  Below,  I will  use  the  superscripts  R and  2 in  order  to 
distinguish  formulas  referring  to  Rouse  and  Zimm  dynamics. 

2.2  Configuration  of  Polymers  at  the  Membrane  • Channels 

The  concentration  of  polymers  directly  at  the  membrane  is  an  important 
parameter  in  my  problem,  it  depends  on  bulk  concentrations  (p0  and  0)  on  both 
sides  of  the  membrane  and  on  the  chemical  potential  of  the  membrane.  It  is 


n 

shown  in  Appendix  A that  this  concentration  is  exactly  one  half  of  the  equilibrium 
concentration  pe„  which  would  be  established  on  the  membrane  when  the  bulk 
concentration  on  both  sides  of  the  membrane  is  p0. 

Pm  = ^(p».)«,  (2.1) 

The  polymers  at  the  membrane  have  different  mobility  depending  on  the  config- 
urations of  these  polymers.  I will  consider  the  How  IT  through  the  membrane  as 
a sum  over  contributions  from  different  channels,  IT  = £ II*.  Each  channel  with 
index  k describes  the  motion  of  polymers  which  cross  the  membrane  k times. 
The  density  of  polymers  in  the  i-th  channel  is  pt  and  Ept  = pm.  The  poly- 
mer density  p,  in  each  channel  can  be  related  to  the  density  in  the  equilibrium 
problem: 

Pi  = (2.2) 

The  surface  density  of  these  polymers  (density  per  unit  area  of  the  membrane) 
can  be  obtained  by  multiplying  pt  by  the  thickness  of  the  layer  of  polymers 
Rs  = b i/w  (see  Fig.  2.1  a).  In  every  time  interval  (f  (or  tf)  the  polymers  in  the 
channel  k are  renewed:  some  polymers  leave  the  channel,  other  polymers  move  in. 
The  contribution  of  each  channel  to  the  transport  per  unit  area  of  the  membrane 

= (2.3) 

where  t,  is  either  the  Rouse  (if)  or  Zimm  ((f)  relaxation  time,  n*  has  the 
dimension  of  [number  of  polymers/(unit  area  x sec)].  The  total  transport  through 


the  membrane  measured  in  (p0Du/b)  units 


n*  - ftf-bS  < > n*  _ jv-o.5  < 1fk2  > (2.4) 

where  < l/k  >,  < l/k2  > are  properly  averaged  negative  moments  of  k : 

<,A>-Eoi.  <ve>-Egi  M 

(The  mass  transfer  formula  can  be  obtained  by  multiplying  nt  by  the  index  of 
the  polymerisation  N.) 

2.3  Calculation  of  Negative  Moments  of  k 

The  average  in  Eq.  2.5  can  be  done  exactly  with  the  use  of  the  equilibrium 
density  distribution  (sec  Eqs.  2.1  and  2.2).  At  any  particular  value  of  the  mem- 
brane potential  p,  the  partition  function  of  the  chain  is  the  sum  of  the  partition 
functions  with  different  numbers  of  crossings  k. 

Z = = (2.6) 

Each  particular  Zk  can  be  factorized  into  the  product  of  two  terms:  Z%  and 
e where  is  the  partition  function  at  p — 0 and  the  relative  weight  e-'** 
contains  all  dependence  on  the  chemical  potential.  The  sum  in  Eq.  2.5  which 
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gives  < 1/1:  > can  be  obtained  by  integrating  Eq.  2.6 

< V* >=  Y°  \ kZ*e  '*  = T°i  (2.7) 

The  sum  in  Eq.  2.5  which  gives  < 1 /*»  > can  be  obtained  by  double  integration: 

< V*J  >=  h WW'dn'  (2.8) 

The  detailed  calculation  of  the  partition  function  Z(p,  N)  is  given  in  Appendix 
B.  This  calculation  is  made  in  a way  similar  to  the  calculation  of  the  partition 
function  in  the  adsorption  problem  [20,21],  I find 

ZfoN)  = t-“  e-^lErM^SN/i))'  (2.9) 

with  Z‘(N)  = 1.  Equations  2.5,  2.7-2.0  give  a closed  form  description  of  the 
flow  through  the  membrane.  The  results  for  various  limiting  cases  can  be  ob- 
tained knowing  that  the  partition  function  Z(p,Af)  is  a universal  function  or 
C = y/iSN/ 4 with  an  additional  factor  e'".  This  factor  means  that  any  config- 
uration has  at  least  one  crossing  of  the  membrane.  For  ( » 1,  Z{i i,  N)  has  an 
asymptotic  form  ~ p~!/Ve“".  Integrating  this  function  I obtain 

< 1/fc  >~  i,  < 1 /**  >~  log(p)  (2.10) 

(and  < 1/*  >~  y/N,  < l/e  >-  log(A')  at  small  p < l/x/77.) 
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2.4  Qualitative  Discussion 

2.4.1  Neutral  Membrane  [n  = 0). 

The  configuration  which  crosses  the  membrane  only  once  (shown  on  the  Fig. 

2.1  b)  has  small  statistical  weight:  each  half  of  the  chain  to  the  right  and  left  of 
the  membrane  can  be  thought  of  as  a «-  ;V/2  step  random  walk,  which  starts  on 
the  membrane  and  does  not  cross  the  membrane  again.  The  probability  of  such 
a random  walk  is  of  order  1/y/N  (see  Eq.  B.4  in  Appendix  B).  The  statistical 
weight  of  a configuration  with  two  such  random  walks  is  ~ 1/N. 

The  average  number  of  crossings  k in  the  typical  configuration  of  a polymer 
chain  is  of  order  of  y/ff.  It  is  shown  in  Appendix  B,  that  the  relative  weight 
of  configurations  with  k crossings  is  — k/N.  In  the  case  of  Rouse  dynamics  it 
happens  that  all  >/N  channels  give  the  same  contribution  to  the  total  transport: 
n£  = W-a,s,  and  the  total  transport  is  (*,„„  = y/N): 

n"  = £ n»  = AT*  (2.11) 

For  Zimm  dynamics  channels  with  small  k give  contribution  proportional  to  l/k: 
nf  = k~lN~l*.  Integration  over  all  k gives  the  logarithmic  N dependence  of  n. 

nz  = logiv 


(2.12) 


2.4.2  Repulsive  Membrane  (p  > 0). 


In  this  case  there  is  a factor  e in  each  11*.  If  /*  < 1/y/N,  this  factor  does 
not  affect  configuration  statistics.  If  p > 1/n/jV,  it  will  produce  an  effective  cut 
off  in  the  summation  over  k in  Eq.  2.10.  The  maximal  number  of  crossings  k is 
of  order  of  k ~ 1/p.  At  p > 1 only  configurations  with  one  crossing  (see  Fig.  2.1 
b)  are  relevant.  For  the  Rouse  model,  the  main  effect  of  membrane  potential  is 
in  controlling  the  number  of  channels  contributing  to  the  transport.  For  Zimin 
dynamics  the  main  contribution  comes  from  small  k.  Instead  of  Eqs,  2.11  and 
2.12  we  have 

n"  = jr'tr**,  n*=  Ar’-'iog^)  (2.13) 

At  p ~ 1 there  is  a crossover  in  both  formulas  to  the  single  channel  regime. 

n«  _ jy-M  n*  = AT1-5  (2.14) 


2.4.3  Attractive  Membrane  (p  < 0). 

At  |p|  > 1/ZN  a polymer  is  adsorbed  on  the  membrane.  The  typical  poly- 
mer configuration  can  be  considered  as  a chain  of  “adsorption  blobs,"  each  blob 
consisting  of  Ng  = l/pa  monomers.  On  the  scale  less  than  the  blob  size  the  chain 
has  the  same  statistics  as  a random  walk  and  crosses  the  membrane  kg  = l/|p| 
times,  so  that  the  adsorption  energy  of  a blob  is  of  order  one  (in  kaT  units). 
The  average  number  of  membrane  crossings  is  proportional  to  the  length  of  a 


polymer: 


k = (N/N,)k,  = |p|iV  (2.15) 

The  factor  e*1''1  increases  the  relative  weight  of  configurations  with  large  k.  The 
main  contribution  to  the  transport  comes  from  about  \/\p\  channels  with  k £ k: 


(2.16) 


In  the  case  of  Zimin  dynamics  there  are  two  main  contributions  to  the  transport, 
one  from  small  k (dis-adsorbed)  configurations  and  the  other  from  large  k (ad- 
sorbed) configurations  at  higher  p.  The  total  transport  is  the  sum  of  these  two 
contribution  n*  = Ilf  + nf 

flf  = * wJ°g>X  [small  *).  Ilf  = [large  k ) (2.17) 

When  e"I''/(|ft|iVJ'a)  ~ 1 (R),  ef‘’K/(ii,Nu)  ~ 1 (Z) , the  overall  permeabil- 
ity of  the  membrane  becomes  comparable  with  the  permeability  of  the  solvent: 
Af-1*5.  In  this  regime,  the  membrane  is  transparent.  This  effect  can  only  be 
observed  in  the  limit  of  very  low  bulk  concentration  of  polymers  pa.  The  concen- 
tration of  polymers  on  the  membrane  is  times  larger  than  pa  and  should  still 
be  less  than  the  concentration  at  which  interactions  between  polymers  become 
important.  Otherwise  the  effects  of  interaction  of  polymers  will  control  the  access 
to  the  membrane  and  the  flow  through  the  membrane  will  be  saturated  or  even 
clogged. 


The  effects  of  excluded  volume  interactions  are  not  considered  in  our  discus- 
sion, but  in  mean  field  approximation  these  effects  can  be  easily  incorporated  in 
our  description,  by  assuming  that  there  are  two  different  chemical  potentials  on 
both  sides  of  the  membrane,  p\  > (12-  The  difference  of  chemical  potentials  is 
proportional  to  the  difference  in  monomer  concentrations  piN  and  P2N  on  both 
sides  of  the  membrane: 


Px  - = (/>,  - p,)NB  = SpNB  (2.18) 

where  B is  the  second  virial  coefficient. 

In  this  case  there  is  a difference  of  the  osmotic  pressure  on  both  sides  of  the 
membrane.  This  results  in  a force  on  the  monomers  when  they  cross  the  mem- 
brane. The  force  is  directed  towards  the  dilute  phase  and  produces  an  average 
drag  velocity 

VR  = D0SpB,  Vz  = DJpNB/k  (2.19) 

If  the  time  tan/t  a polymer  needs  to  cross  the  membrane  is  shorter  than  the 
diffusion  time  tR(tR),  the  transport  If  = p/tRg  ~ p1  is  nonlinear.  This  happens 
under  the  condition  that 

N'*{pbB)  < 1 (2.20) 

The  configurations  with  multiple  crossings  of  the  membrane  are  suppressed 
in  this  regime  because  of  the  low  probability  of  returns  to  the  dense  phase. 
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2.5  Discussion 

I considered  a problem  of  a polymer  chain  penetrating  through  a porous  mem- 
brane by  calculating  probabilities  of  different  configurations  of  a polymer  and  the 
corresponding  mobilities.  I found  the  optimal  configurations  which  give  the  main 
contrioution  to  the  transport  through  the  membrane.  In  case  of  Rouse  dynamics, 
when  p = 0 the  configurations  with  different  numbers  of  membrane  crossings  k 
give  equal  contribution.  At  nonzero  p these  contributions  arc  proportional  to 
e~uk,  which  give  preference  to  configurations  with  small  (large)  k for  a repulsive 
(attractive)  membrane. 

If  hydrodynamic  interactions  are  important,  the  main  friction  effects  during 
the  motion  of  a polymer  are  localized  at  the  points  where  the  polymer  crosses 
the  membrane.  In  this  case,  the  main  contribution  to  the  transport  comes  from 
configurations  with  small  numbers  of  k.  If  the  membrane  is  attractive,  the  main 
contribution  is  switched  to  configurations  with  large  k. 

The  process  of  diffusion  through  the  membrane  can  be  used  for  separation  of 
polymers  with  different  lengths.  The  effectiveness  of  this  process  can  be  character- 
ized by  the  dependence  in  the  expressions  for  the  transport  Eqs.  2.11-2.14. 
The  exponent  a varies  from  0.5  to  2.5.  There  is  also  an  additional  exponential 
dependence  on  iV  in  Eqs.  2.15  and2.16  for  the  attractive  membrane.  The  use 
of  adsorbing  membranes  (p  < 0)  might  not  be  very  practical  because  of  low 
permeability  and  the  possibility  of  changing  membrane  properties  by  adsorbed 
polymers. 


CHAPTER  3 


DETERMINISTIC  MODEL  OF  DNA  GEL  ELECTROPHORESIS  IN  A 
STRONG  ELECTRIC  FIELD 

3.1  Beyond  the  Tube  Model 

The  motion  of  a high  molecular  weight  polymer  in  a gel  is  constrained  by  the 
strands  of  the  gel.  The  strands  limit  the  displacement  of  a polymer  in  the  direc- 
tion transverse  to  its  contour.  In  the  tube  model,  these  topological  constraints 
of  network  strands  arc  replaced  by  a confining  tube  (18, 19]  (of  diameter  similar 
to  the  size  of  the  pores).  In  the  absence  of  an  external  field,  the  chain  moves 
through  a gel  along  this  tube  by  a reptation  mechanism  [1,2].  The  presence  of  an 
electric  field  introduces  a bias  in  the  motion  of  charged  polymers.  The  thermally 
activated  motion  along  the  tube  in  the  field  direction  is  more  preferable  than  in 
the  opposite  direction.  The  reptation  with  a bias  is  the  dominant  mechanism 
of  polymer  motion  in  relatively  weak  fields  [22-30].  This  model  works  remark- 
ably well  for  the  intermittent  regimes  even  at  large  fields  [27].  As  the  electric 
field  E becomes  stronger,  the  chain  can  no  longer  move  exclusively  along  its  lin- 
ear contour  (tube).  If  the  electric  energy  of  the  charged  entanglement  strand  is 
comparable  to  the  thermal  energy  k0T 


qEa  as  kgT 


(3-1) 


the  walls  of  the  tube  are  no  longer  able  to  confine  the  polymer.  Here  q is  the 
charge  per  “blob.”  (The  section  of  the  chain  contained  in  one  pore  of  size  a is 
called  a “blob.”) 

The  external  field  can  pull  double  folded  loops  of  the  molecule  through  the 
gaps  between  the  gel  fibers  (31).  The  strands  of  the  chain  leak  out  of  the  tube 
in  the  shape  of  hernias  under  the  influence  of  the  field.  These  hernias  have 
been  found  both  in  computer  simulation  [31-35],  and  by  direct  observation  of 
individual,  fiuorescently-stained  DNA  molecules  under  the  microscope  [36-39]. 
There  is  an  indication  that  the  average  number  of  hernias  increases  with  molecular 

A remarkable  tree-like  structure  of  a moving  polymer  was  found  in  the  com- 
puter simulations  by  Duke  and  Viovy  [40, 41]  and  by  Deutscu  [42].  The  results 
of  these  simulations  indicate  that  the  motion  of  an  ultra-high  molecular  weight 
polymer  is  qualitatively  different  from  the  “hang  and  slip”  mechanism  of  motion 
of  a lower  molecular  weight  polymer. 

In  this  chapter,  I present  results  for  the  recently  proposed  "deterministic 
model"  [43]  of  DNA  motion  through  a gel  in  a strong  electric  field.  I obtain  de- 
tailed information  about  the  statistical  properties  of  loops  and  find  novel  thresh- 
old features  in  their  dynamics.  The  model  is  almost  deterministic  and  allows  the 
study  of  DNA  electrophoresis  for  molecular  weights  unreachable  in  usual  simula- 
tion. In  the  next  section,  I describe  the  deterministic  model  of  gel  electrophoresis. 
In  section  3.3,  the  analytical  predictions  for  this  model  will  be  presented.  In  sec- 
tion 3.4,  I give  the  results  of  my  computer  simulation.  In  the  last  section  of 


this  chapter,  I 


the  limitations  of  the 
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3.2  Deterministic  Model 

A schematic  sketch  of  a typical  configuration  of  a moving  linear  polymer  with 
many  loops  (hernias)  is  shown  in  Fig.  3.1.  It  is  an  oversimplified  sketch,  and 
the  displacements  perpendicular  to  the  electric  field  are  considered  unimportant 
and  arc  not  shown  to  scale.  The  typical  configuration  consists  of  a hierarchy  of 
hernias  of  various  sizes.  The  polymer  is  fully  stretched  everywhere  except  near 
the  tips  of  the  hernias,  where  the  linear  tension  is  zero  (35).  I assume  that  at  the 
tips  of  the  hernias  the  DNA  is  in  the  coil-like  configuration,  with  a small  amount 
of  stored  length  that  is  negligible  in  my  simulations.  Hernias  engage  in  a tug  of 
war  with  each  other  and  those  that  are  down  the  field  grow  at  the  expense  of  the 
others  (44).  While  these  hernias  grow,  their  tips  encounter  strands  of  the  gel  and 
might  split,  creating  new  hernias.  Some  of  these  small  hernias  grow,  while  others 
disappear. 

The  motion  of  the  large  sections  of  the  chain  (for  example  section  AOB  of 
Fig.  3.1)  can  be  described  by  simple  deterministic  equations.  The  tension  at 
the  hernia  tips  is  zero.  This  implies  that  instantaneously  there  is  no  mechanical 
force  between  the  strand  and  the  rest  of  the  polymer.  This  section  of  the  chain 
is  moving  because  of  the  asymmetry  of  the  lengths  AO  and  OB.  The  equation 


of  motion  for  this  section  is 


(3-2) 


where  { is  the  friction  coefficient  per  chain  section  of  size  a and  with  charge  q, 
v0  = qE/£  and  Lao  and  Lob  are  the  lengths  of  sections  >10  and  OB,  respectively. 
The  velocity  of  this  segment  depends  only  on  the  asymmetry  factor  AAOb  = 
(Lao  — Lob) /(Lao  + Lon)  of  the  two  competing  strands  .40  and  OB  hanging 
on  the  common  gel  strand  "hook”  O.  If  one  knows  the  velocities  of  each  strand, 
it  is  easy  to  write  the  velocity  of  the  tip  points: 


\(aaob  Abba),  = "o  - Aaob)-  (3.3) 


Because  of  the  strong  field,  the  chain  is  fully  stretched  most  of  the  time  and 
the  total  length  of  the  polymer  (the  sum  of  length  of  all  strands)  remains  constant. 
I extended  this  description  to  the  cases  when  some  parts  of  the  polymer  are  in 
a free  fall.  This  happens  when  one  of  the  ends  (on  the  right  or  left  side  of  Fig. 
3.1)  is  pulled  over  the  (entanglement)  hook  by  the  rest  of  the  chain.  1 make  the 
assumption  that,  when  this  happens,  the  strand  moves  with  a homogeneous  free 
fall  velocity  e„. 

Equations.  3.2  and  3.3,  written  for  every  strand  of  a polymer  chain,  form 
a set  of  coupled  linear  differential  equations.  These  equations  can  be  solved 
numerically.  When  the  length  L,  of  a strand  becomes  zero  this  strand  disappears 
and  the  set  of  equations  of  motion  for  the  adjacent  strands  is  rewritten.  The  above 


Figure  3.2:  The  schematic  picture  of  the  tip  splitting.  After  traveling  the  pore 
sac  distance  a,  the  tip  could  split  into  two  with  probability  p.  During  this  process 
the  total  length  of  the  chain  does  not  change. 

equations  are  completely  deterministic,  but  they  describe  only  the  evolution  of 
existing  hernias  and  their  disappearance.  The  creation  of  new  hernias  must  be 
added  to  this  model. 

I assume  that  each  time  a tip  travels  distance  a and  moves  from  one  pore  to 
the  neighboring  one,  there  is  a probability  p that  it  will  hit  an  obstacle  and  will 
split  into  two  (see  Fig.  3.2).  Two  small  hernias  with  length  I,,  and  lj  are  created 
as  a result  of  the  split. 

The  total  length  is  conserved  during  this  process  (see  Fig.  3.2): 


(3.4) 
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The  initial  sizes  of  new  hernias  (li.li)  are  determined  randomly  in  the  interval 
(0,  ^j)  in  units  of  the  pore  size  a.  This  is  the  only  uondeterministic  part  of  my 
model.  After  the  creation  of  the  new  hernias,  their  dynamics  is  controlled  by 
Eqs.  3.2  and  3.3. 

I study  the  dynamics  of  linear  and  circular  polymers  with  length  ranging 
between  102  and  105a,  which  corresponds  to  10_1  - 102  megabases.  Everywhere 
below  I will  use  pore  size  units  (a  typical  pore  size  is  1000  A and  containing 
approximately  1000  base  units,  see  for  example  Duke  and  Viovy  [41]).  I compare 
the  results  for  linear  and  circular  polymers  in  order  to  isolate  and  to  understand 
the  end  effects. 

The  total  number  of  linear  differential  equations  to  be  integrated  for  each  time 
step  is  twice  the  number  of  hernias.  This  number  grows  only  logarithmically  with 
molecular  weight.  The  evolution  time  of  a DNA  configuration  is  the  typical  time 
during  which  the  polymer  is  displaced  along  the  field  by  its  average  size.  This 
time  is  proportional  to  the  molecular  weight  of  the  polymer.  Thus,  the  actual 
computation  time  of  my  simulation  grows  only  as  .V  log  N with  the  molecular 
weight  N.  The  simulation  for  the  1 megabase  DNA  (overall  length  1000  pore 
sizes)  for  about  100  evolution  times  takes  a few  minutes  on  SUN  Sparc  690  (33 
MHz  CPU  speed)  workstations. 

First,  I consider  a circular  polymer.  For  a linear  polymer  the  results  arc 
basically  the  same,  but  with  broader  statistical  distributions. 


3,3  Analytical  Results 


In  Obukhov  and  Rubinstein  (43),  It  is  argued  that  the  average  drift  velocity 
is  independent  of  the  length  of  the  DNA.  It  is  also  explained  the  origin  of  the 
statistical  self-similarity  of  the  moving  polymer.  Both  of  these  properties  arc 
due  to  the  [act  that  the  main  equations  controlling  the  evolution  of  a particular 
configuration  are  scale  invariant.  The  velocity  of  a particular  strand  calculated 
with  Eq.  3.2  depends  only  on  the  ratio  of  the  geometrical  sixes  LAo,  Lob , but 
not  on  their  absolute  values.  I can  study  the  dynamics  of  some  configuration  of 
a polymer  of  length  N using  Eqs.  3.2  .aid  3.3,  and  then  multiply  all  longitudinal 
hernia  sizes  and  times  by  an  arbitrary  factor  A.  I obtain  the  description  of  the 
dynamics  for  the  polymer  of  size  A1V.  The  only  process  which  violates  this  scaling 
is  the  microscopic  mechanism  of  the  creation  of  new  hernias,  and  I will  address 
this  issue  in  the  numerical  simulation  section.  The  lifetime  of  each  side  branch 
of  a self-similar  tree-like  structure  proposed  by  Obukhov  and  Rubinstein  [43]  is 
proportional  to  its  distance  from  the  leading  tip  (Fig.  3.1).  This  is  the  time  during 
which  this  branch  (that  was  created  in  the  vicinity  of  a leading  tip)  competed 
with  others  and  survived,  while  the  tip  moves  forward.  When  all  velocities  in  the 
problem  are  proportional  to  the  size  of  the  side  branch  should  be  of  the  same 
order  of  magnitude  as  the  distance  from  the  leading  tip.  I expect  that  on  each 
size  scale  the  number  of  surviving  branches  will  be  the  same  (of  order  unity); 
hence,  the  total  number  of  hernias  n is  proportional  to  the  logarithm  of  the  total 


length  ;V  of  the  DNA. 


n ~ logJV.  (3.5) 

The  probability  to  End  a hernia  in  the  size  interval  (L,  L + SL)  is 

P(L)~1/L.  (3.6) 


3.3.2  Critical  Velocity 

The  dynamics  of  the  tree-like  structure  (Fig.  3.1)  is  delicately  related  to  the 
process  of  creation  of  new  branches  (hernias)  and  their  disappearance.  When  the 
tip  of  a hernia  splits  in  two  (see  Fig.  3.2),  the  new  baby  hernias  are  trying  to 
outgrow  each' other.  1 assume  that  the  velocities  of  both  hernias  0i  and  Vj  depend 
on  the  large  scale  details  of  the  configuration  and  are  constant  during  this  process. 
We  also  assume  that  the  asymmetry  parameter  A = (li — h)/(li +fz)  of  the  initial 
sizes  l|,/j  of  new  hernias  is  uniformly  distributed  in  the  interval  (-1,1).  If  the 
velocity  of  the  splitting  tip  is  negligible  or  zero  (»i  sz  vs  c:  0),  the  time  needed  for 
the  larger  of  the  two  hernia  to  overweigh  the  other  is  very  short  and  is  of  order 
Is)/!),.  In  the  other  limit,  when  tire  tip  is  moving  with  free  fall  velocity 
U|  = Uj  = n„,  and  splits  into  two,  the  lifetime  of  the  new  hernias  is  infinite.  This 
is  because  the  flow  of  mass  between  the  two  baby  hernias  in  Eq.  3.1  v is  always 
smaller  than  the  free  fall  velocity,  when  the  input  of  mass  to  each  hernia  is  o„. 


In  Appendix  C,  I demonstrate  that  there  exists  a special  ratio  of  the  sizes  of 
baby  hernias 

A m ‘‘-k  = ("■  ~ (3  7) 

Ii  + (a  (2«o-(vi+»a)) 

at  which  in  the  flow  of  mass  between  the  two  hernias  is  balanced  in  such  a way, 
that  the  asymmetry  ratio  A^ag  remains  constant.  In  the  vicinity  of  Aging  the 
lifetime  t is  a diverging  function  of  1-4  — 

In  the  Appendix  C,  I calculate  the  lifetime  t of  a new  hernia,  averaged  over  all 
possible  asymmetries  -1  < A < 1.  The  remarkable  new  result  of  this  calculation 
is  that  there  is  a special  critical  velocity  of  the  tip, 

vc  = v„/2  (3.8) 

at  which  the  average  lifetime  diverges.  If  the  velocity  of  the  tip  is  less  then  ve, 
the  average  lifetime  of  new  hernias  is  finite,  if  v > ve  then  this  time  is  infinite.  It 
means  that  statistically  at  velocities  v > ve  branching  will  lead  to  formation  of 
stable  new  hernias. 

3.3,3  Critical  Ratio 

Suppose  that  there  is  a simple  configuration  of  a large  circular  DNA  (Fig. 
3.3).  The  right  part  slowly  overweighs  the  left  part  and  small  hernias  are  created 
from  time  to  time  at  the  right  tip  of  this  structure.  These  small  hernias  have 
symmetric  input  (i.'i  = vg)  from  the  main  strand.  Their  longevity  is  controlled 
by  the  velocity  of  the  tip  v — [vi  + vg)/2.  If  v < vc  they  will  disappear,  if  u > ve, 
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Figure  3.3:  a),  b)  Critical  ratio  for  a V-shape  configuration  for  a circular  polymer. 
One  branch  is  overweighing  the  other.  When  the  ratio  of  the  length  of  two 
branches  Lj  : L\  exceeds  3:1,  the  stable  branching  of  the  tip  becomes  possible, 
c),  d)  For  a linear  polymer,  the  critical  ratio  is  ij  : L\  = 5:1. 
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the  small  hernias  will  start  to  grow.  The  velocity  v depends  on  the  ratio  of  sizes 
of  the  right  and  left  parts  of  the  DNA,  Li  and  Lj.  Using  Eq.  3.2 1 find  that  there 


at  which  the  velocity  of  the  moving  tip  is  critical.  If  this  ratio  is  smaller  than 
3:1,  then  v < ve  and  all  small  hernias  will  disappear.  As  one  part  overweighs 
the  other,  the  ratio  Lf.L,  grows  with  time.  When  the  ratio  becomes  larger  than 
3:1,  the  velocity  of  the  tip  exceeds  the  critical  velocity  » > vc.  Small  hernias, 
created  at  the  tip,  become  stable  and  start  to  grow.  This  critical  ratio  was  clearly 
observed  in  my  simulation  of  circular  polymers. 

3.3.4  Linear  Polymer 

The  calculation  scheme  presented  in  Appendix  C can  be  easy  modified  for  the 
case  of  a branching  of  a moving  end  of  a linear  polymer;  see  Fig.  3.3  (c)  and  (d). 
In  this  case  I obtained  for  the  critical  velocity  of  a leading  tip: 


The  critical  asymmetry  ratio  at  which  the  branching  of  a leading  tip  will  occur 


: l,  = 5 : 1 


(3.11) 
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Figure  3.4:  Average  number  of  hernias  as  a function  of  molecular  site  of  the 
polymer.  It  shows  logarithmic  molecular  size  dependence  for  both  circular  and 
linear  polymers  (for  N > 1000  ). 


3.4.1  Average  Number  of  Hernias 

1 have  performed  simulations  in  the  size  N range  from  102  to  10s  units  (this 
corresponds  to  10*  ~ 10s  base  pairs)  for  circular  and  linear  polymers.  I found 
very  slow  dependence  of  the  average  number  of  hernias  on  molecular  weight  (see 
Fig.  3.4). 

For  large  (N  > 1000)  I observed  a logarithmic  behavior  of  the  average  number 
of  hernias  (see  Fig.  3.4)  for  both  linear  and  circular  polymers.  This  confirms  the 
theoretical  prediction  [43]  of  the  statistical  self  similarity  of  the  moviug  polymer 
chain.  This  result  is  in  qualitative  agreement  with  the  results  of  Duke  and  Viovy 
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Figure  3.5:  Average  velocity  of  center  of  mass  as  a function  of  molecular  size 

N.  The  velocity  is  independent  of  the  molecular  size  of  polymer  up  to  possible 
logarithmic  corrections.  The  average  velocity  of  a linear  chain  is  larger  than  that 
of  a circular  chain  due  to  the  end  effect. 

[40),  who  observed  ramified  structures  in  simulation  of  megabase  long  chains 
(N  > 1000).  For  shorter  chains  (IV  < 1000),  the  statistical  self-similarity  breaks 
down,  and  I do  not  observe  tree-like  structures.  The  polymer  spends  most  of  its 
time  in  a simple  configuration  with  only  two  branches. 

3.4.2  Average  Velocity 

The  results  for  the  average  velocity  of  a moving  polymer  at  long  times  are 
shown  in  Fig.  3.5.  For  the  circular  polymer  the  average  velocity  approaches 

O. 31c,,.  For  a linear  polymer  the  average  velocity  is  0.34c,,. 

The  result  that  the  average  polymer  velocity  is  independent  of  the  molecular 
weight  is  consistent  with  the  model  [43]  and  numerical  simulation  [40].  The  very 
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slow  decrease  in  the  velocity  cau  be  explained  by  the  logarithmic  corrections 
~ 1 / log  ;V  to  the  average  velocity  because  the  velocity  is  averaged  over  both 
non-tree-like  configurations  and  tree-like  configurations. 

3.4.3  Periodicity  of  the  Motion  of  a Circular  Polymer 

I find  strong  periodicity  in  the  motion  of  circular  DNA  through  the  gel.  The 
number  of  hernias  and  velocity  of  the  center  of  mass  plotted  as  functions  of  time 
(Fig,  3.6  upper  and  lower  graph  respectively)  show  this  periodic  pattern.  The 
average  time  of  a cycle  increases  linearly  with  the  molecular  weight  of  DNA  (Fig. 
3.7). 

Successive  interchanges  between  tree-like  structure  and  V-shape  configuration 
take  place  regardless  of  molecule:  size.  V-shape  configurations  can  be  recognized 
as  small  plateaus  with  n = 2 on  the  upper  graph  and  segments  of  exponentially 
growing  velocity  on  the  lower  graph  of  Fig.  3.6. 

The  series  of  pictures  in  Fig.3.8  shows  a typical  evolution  pattern  within  a Cy- 
cle. The  cycle  starts  with  a V-shape  configuration  with  some  asymmetry  between 
left  and  right  branches.  This  asymmetry  increases  slowly  until  the  asymmetry 
ratio  exceeds  3:1.  At  this  moment  the  branching  starts  at  the  leading  tip.  The 
further  evolution  can  be  described  in  terms  of  a statistically  self-similar  tree-like 
structure  due  to  the  competition  between  branches  and  the  creation  of  new  ones. 
If  the  two  largest  branches  of  the  tree  have  comparable  length  (Fig.3.8(d-f)), 
the  branching  in  each  of  them  slows  down  and  smaller  branches  start  to  disap- 
pear. If,  after  the  disappearance  of  smaller  branches,  the  length  ratio  of  two  large 
branches  remains  less  than  3:1,  the  situation  returns  to  the  beginning  of  the  cycle 
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Figure  3.6:  Number  of  hernias  and  center  of  mass  velocity  as  a function  of  time 
for  a circular  polymer.  The  figure  shows  a periodic  pattern.  A large  number  of 
hernias  are  indicative  of  tree-like  structures.  The  plateaus  on  the  upper  graph 
with  n = 2 indicate  a simple  V-shape  configurations.  Only  the  hernias  larger 
than  5a  are  counted.  The  velocity  of  the  center  of  mass  is  0.25tio  at  the  onset 
of  the  growth  in  hernia  numbers  (this  corresponds  to  the  critical  linear  velocity 
Vc  = 0.5v„). 


Figure  3.7:  Dependence  of  the  period  on  the  molecular  weight  of  the  circular 
polymer.  One  period  is  defined  as  the  time  between  two  consecutive  major  peaks 
in  the  velocity  of  the  center  of  mass  (see  Fig.  3.6  and  Fig.  3.9).  The  data  ore 
averaged  over  100  periods. 

(Fig.3.8(f)).  The  duration  of  each  stage  of  this  process  is  about  the  same  and 
is  proportional  to  the  total  length  of  the  DNA,  During  the  simulation,  it  is  only 
possible  to  clearly  identify  stage  I,  (plateau  on  Fig.3.6).  The  boundary  between 
stages  II  and  III  (increasing  and  decreasing  velocity  of  center  of  mass)  is  often 
ill-defined. 

The  circular  DNA  stays  in  a simple  V-shape  configuration  about  30%  of  total 
time  (see  Fig.3.6).  I observe  similar  behavior  in  the  case  of  a linear  chain.  But 
the  cycle  for  linear  chains  is  not  as  regular  as  for  circular  chains.  I explain  this 
difference  by  the  “free  ends  effect1' (see  Fig  3.1).  When  the  chain  is  in  the  tree-like 
regime  almost  all  of  its  branches  are  double  folded.  When  two  large  branches  of 
comparable  length  are  formed,  they  are  also  often  double  folded.  The  falling 
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Figure  3.8:  Time  evolution  of  a circular  polymer.  The  figures  show  a typical 
cycle  between  two  consecutive  V-shape  configurations.  The  leading  tip  starts  to 
branch  as  soon  as  the  ratio  of  lengths  of  the  two  branches  exceeds  the  critical 
ratio  3:1,  (b).  The  tree  is  formed  (c-d)  and  then,  after  formation  of  two  major 
competing  branches,  is  simplified  again  (e-f). 
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Figure  3.9:  Number  of  hernias  and  center  of  mass  velocity  as  a function  of  time 
for  a linear  polymer.  The  periodicity  is  much  less  clear,  then  on  Fig.  3.6  due  to 
free  failing  ends. 

free  ends  may  create  an  asymmetry  of  the  configuration,  which  is  larger  than  the 
critical  one,  and  the  leading  hernia  starts  branching  again  . I found  that  due  to 
the  end  effect  a linear  polymer  stays  in  a simple  V-shape  configuration  only  10% 
of  the  time  (Fig.3.9). 


3.4.4  Tip  Splitting  Probability  and  Size  Distribution 

Most  of  my  computer  simulation  were  done  with  the  probability  of  tip  splitting 
p = 0.5  (each  time  a tip  travels  the  distance  of  the  pore  site  it  has  a 50%  chance  of 
splitting).  1 checked  the  sensitivity  of  my  model  to  the  choice  of  this  parameter. 
In  Fig.3.10  the  hernia  size  distribution  function  P(L)  is  plotted  as  a function  of 
L on  logarithmic  scale  for  different  splitting  probabilities  p ranging  from  0.005 
to  1.  It  clearly  demonstrates  the  universality  of  the  law  P{L)  ~ l/L  (see  Eq. 
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Figure  3.10:  Size  distributions  of  hernias  plotted  for  a Megabase  (N  = 1000) 
circular  polymer  at  various  tip  splitting  probabilities  p.  The  straight  line  indicates 
l/L  dependence.  At  p > 0.2,  the  curves  are  indistinguishable. 

3.6)  for  p > 0.01.  If  the  probability  of  tip  splitting  becomes  too  low  (p  < 0.01), 
the  distribution  P[L)  is  universal  only  for  small  size  hernias.  The  distribution 
resembles  that  of  short  chains,  where  configurations  with  two  branches  give  the 
dominant  contribution  to  the  statistics  of  large  hernias  with  L ~ rV. 

On  Fig.3.11  the  frequency  of  hernias  of  unit  size  P(l)  is  shown  for  different 
tip  splitting  probabilities  p ranging  from  0.005  to  1.  For  the  improvement  of 
the  accuracy,  the  frequency  P(l)  is  calculated  as  the  intercept  of  curves  shown 
in  Fig.3.11  with  the  line  L = 1.  The  slope  1/2  in  the  double  logarithmic  plot 
suggests  a following  scaling  behavior  at  constant  L\ 

PW~p'n.  (3.12) 
The  frequency  P saturates  and  becomes  constant  at  p > 0.2,  1 have  no  analytical 


Figure  3.11:  Values  of  Che  hernia  distribution  function  for  unit  size  hernias  P(L  = 
1)  as  a function  of  tip  creation  probability  p on  logarithmic  scale.  The  slope  of 
this  plot  suggests  a power  law  P(L  = 1)  ~ pV2. 
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explanation  for  the  observed  exponent  1/2.  I can  combine  the  scaling  laws  Eq. 
3.6  and  Eq.  3.12  in  a single  phenomenological  formula: 

P{L)~p'l’/L  /(JV/p1/J)  (3.13) 

where  f(x)  is  the  cross-over  function,  /(x)  = x for  small  x and  /(x)  = const  for 


3.5  Summary  and  Discussion 

I presented  new  analytical  and  numerical  results  for  a simple  deterministic 
model  of  inegabase  DNA  undergoing  in  strong  gel  electrophoresis  field.  The 
results  can  be  summarized  as  follcv.s 

(i)  I confirm  the  statistical  self-similarity  of  the  shape  of  a moving  polymer. 
The  average  number  of  hernias  is  proportional  to  the  logarithm  of  the  molecular 
weight  (Eq.  3.5). 

(ii)  The  average  velocity  is  independent  of  molecular  weight  (with  possible 
logarithmic  corrections).  It  approaches  31%  (34%)  of  the  free  falling  velocity  for 
a circular  polymer  (linear  polymer). 

(iii)  1 find  the  periodic  evolution  pattern  similar  to  that  observed  by  Deutsch 
[42)  in  the  motion  of  a circular  polymer  with  a period  proportional  to  its  molec- 
ular weight.  During  a typical  period  the  shape  of  polymer  changes  from  a simple 
V-shape  configuration  to  a growing  tree-like  conformation  and  then  to  a decaying 
tree-like  conformation.  A circular(linear)  polymer  stays  in  the  V-shape  configu- 
ration about  30%  (10%)  of  the  time  (which  is  truly  remarkable).  The  periodicity 
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is  much  less  pronounced  for  a linear  chain  due  to  perturbative  elfect  of  free  ends. 

(iv)  In  my  simulation  I observed  the  theoretically  predicted  threshold  for 
branching  ve  = va/2  and  the  critical  ratio  3:1  for  a circular  polymer. 

(v)  The  statistical  properties  of  a moving  polymer  are  almost  independent  of 
the  probability  p of  tip  splitting  in  the  interval  0.2  < p < 1 because  of  destructive 
interaction  of  numerous  small  hernias  on  small  scales.  At  low  p < 0.2  the  hernia 
size  distribution  scales  with  polymer  size  L and  probability  of  tip  splitting  p (Eq. 
3.13). 

The  main  advantage  of  the  deterministic  model  is  its  simplicity,  which  allows 
fast  simulation  of  ultra-high  molecular  weight  polymer.  A possible  drawback 
of  the  model  is  the  absence  of  stored  length.  The  effects  of  stored  length  are 
negligible  fur  large  polymers  because  of  strong  stretching  during  most  of  the 
cycle.  But  stored  length  can  be  important  during  the  free  fail  regime.  The  effect 
of  elastic  forces  during  the  free  fall  can  change  the  free  fall  velocity  of  different 
parts  of  a polymer  by  a factor  l/N,  which  might  not  be  very  important,  but 
the  accumulation  of  stored  length  during  the  free  fall  can  result  in  a modified 
frequency  of  tip  splitting  [27j.  I can  argue  that  the  increased  frequency  of  tip 
splitting  might  not  be  noticeable  because  of  the  saturation  effect,  but  this  can 
only  be  checked  by  a future  direct  simulation  of  a more  complicated  model. 


CHAPTER  4 


CHARGED  POLYMERS 

Polyelectrolytes  are  polymers  with  ionizable  groups  [45-51].  When  the  poly- 
mers are  in  polar  solvents,  these  groups  are  dissociated  from  the  polymer  and 
move  away  from  the  polymer.  Therefore,  the  remaining  polymer  become  charged. 
The  dissociated  ion  groups  are  called  “counterions'1.  Polyelectrolyte  systems  are 
very  common  in  nature  and  have  numerous  applications  because  they  are  solu- 
ble in  water.  Water  is  a bad  solvent  for  many  polymers.  But  if  the  molecules 
are  charged,  the  monomers  repel  each  other,  therefore  the  polymer  can  dissolve 
in  the  water.  Most  biopolymers  are  polyelcctrolytes.  Well  known  examples  are 
DNA,  RNA  and  proteins.  A large  uumber  of  synthetic  polyelectrolytes  exist  as 
well.  Polyacrylic  acid  and  polyethylene  sulfide  are  examples.  They  are  used  as 
super-absorber  of  water. 

Polyelectrolyte  solutions  show  a number  of  remarkable  differences  from  that 
of  neutral  polymers  from  the  basic  properties  like  viscosity,  osmotic  pressure  and 
light  scattering  pattern  [52-71). 

Because  of  their  fundamental  relevance  to  biology  and  biochemistry  and  due 
to  their  hydro-solubility,  polyelectrolytes  have  been  of  continuous  interest  in  poly- 
mer science  and  a considerable  theoretical  work  has  been  devoted  to  the  various 
aspects  of  polyelectrolytc  solution  [72-79].  Despite  of  this,  the  understanding  on 
a polyelectrolyte  solution  still  remains  largely  unsatisfactory. 
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The  main  theoretical  difficulty  arises  with  the  long  range  nature  of  the  Coulomb 


interactions.  The  distribution  of  charges  introduces  more  than  one  length  scale, 
which  makes  useful  techniques  for  neutral  polymers  invalid  for  polyelectrolytes. 
The  renormalization  group  theory  or  scaling  theory  are  powerful  techniques  for 
the  neutral  polymer.  In  charged  polymer  systems  it  is  difficult  to  apply  them 
because  of  the  long  range  nature  of  the  Coulomb  potential  [1]. 

Polyelectrolyte  solutions  contain  more  than  one  component.  The  presence  of 

intermediate  length  scale.  The  complicated  interplay  between  small  and  large 
scale  interactions  have  a non  trivial  influence  on  the  local  structure  of  the  poly- 
electrolyte chain.  For  example,  the  local  chain  structure  influences  the  counterion 
condensation,  and  counterion  condensation  in  turn  modifies  the  long  range  inter- 
action. The  long  range  interaction  is  known  to  have  a substantial  contribution 
on  the  stiffness  of  the  chain  [SO]. 

Due  to  these  reasons,  the  charged  polymers  have  been  known  as  an  infamous 
problem  because  of  the  discrepancy  between  theoretical  prediction  and  experi- 
mental results  [80]. 

As  an  example  of  this  controversy,  there  is  a long  standing  debate  regarding 
the  "persistence  length"  [81-84, 88-100, 131],  which  is  the  local  stiffness  of  a 
polymer  chain.  The  persistence  length  is  one  of  the  important  scaling  parameters. 
The  first  scaling  theory  was  provided  by  de  Gennes  et  al  [52].  In  this  model,  it  is 
assumed  that  the  persistence  length  of  the  chain  is  proportional  to  the  screening 
length  which  is  the  interaction  range  of  the  electrostatic  interaction  between 


the  monomers.  In  another  analytic  approach  by  Odijk  and  independently  by 
Skolnick  and  Fixman(OSF)  [81,82],  it  was  predicted  that  the  persistence  length 
is  proportional  to  the  square  of  the  screening  length.  They  calculated  the  energy 
required  to  bend  a charged  polymer  and  found  that  the  stiffness  of  the  chain  is 
strongly  enhanced  by  electrostatic  interactions. 

However,  several  computer  simulations  [94-100]  indicate  that  the  persistence 
length  of  a polyelectrolyte  chain  with  counterions  are  much  shorter  than  the 
OSF  prediction.  The  perturbation  of  ions  surrounding  the  polyelectrolyte  chain 
is  much  stronger  than  the  Debye-Hiickel  estimation.  Therefore,  the  screening  of 
repulsive  interactions  between  charges  on  the  polymers  by  couterions  is  much 
stronger  than  the  screening  given  by  Debye-Hiickel  theory.  In  contrast  to  theory, 
experiments  clearly  indicate  that  the  persistence  length  is  proportional  to  the 
screening  length  (Ip  — 1/i/c ).  As  is  discussed  above,  even  after  much  theoretical 
work  the  origin  of  the  stiffness  of  a charged  polymer  is  not  completely  clear. 

In  spite  of  the  unknown  microscopic  length  scale,  there  has  been  progress  on 
the  study  on  conformational  properties  of  charged  polymers  [5-7. 101, 102]. 

In  a good  solvent,  the  picture  of  a polyelectrolyte  is  a series  of  stretched  elec- 
trostatic blobs  [52]  (see  Fig.  4.1).  Inside  of  each  electrostatic  blob,  the  Coulomb 
repulsion  is  not  sufficient  to  deform  the  chain. 

In  a poor  solvent,  the  chain  collapses  in  the  form  of  a compact  globule  [1, 
103],  The  polymer  chain  has  negative  second  virial  coefficient  which  corresponds 
to  effective  attraction  between  the  monomers.  The  overall  shape  of  a charged 
polymer  in  a poor  solvent  is  determined  by  a balance  of  Coulomb  repulsion  and 
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et  al.  {5]  to  polyelectrolytes  in  a poor  solvent.  I study  the  conformation  of  a 
ring  polyampholytes  in  the  following  chapter.  In  addition  to  the  necklace  like 
deformation  with  respect  to  the  change  of  the  net  charge,  the  ring  polymer  has 


glassy  phase  and  hysteresis  effects. 


CHAPTER  5 


MULTIPLE  CONFORMATIONS  IN  POLYAMPHOLYTES 

The  instability  of  a weekly  charged  polyampholyte  with  respect  to  formation 
of  a necklace  of  globules  was  recently  pointed  out  by  Kantor  and  Kardar  [6,7].  I 
extend  these  results  from  a linear  to  a ring  polymer.  There  is  new  feature  in  a 
ring  polymer  - the  competition  between  linear  and  ring  necklace  formation.  This 
is  analysed  theoretically  using  mean  field  theory  and  numerically  verified.  There 
are  hysteresis  effects  in  the  shape  of  the  clusters.  If  the  total  charge  is  increased, 
the  linearly  extended  configurations  are  dynamically  more  preferable.  When  the 
charge  is  decreased,  the  circular  shapes  are  preferable.  I remark  on  the  probable 
connection  to  the  multiple  phase  transitions  observed  in  polyampholyte  gels. 

5.1  Cascade  of  Transitions  in  a Linear  Polymer 

Polyampholytes  (PA)  are  charged  polymers  which  have  positive  and  negative 
charges  distributed  randomly  along  their  backbone  [110].  The  strength  of  electro- 
static interaction  and  the  relative  ratio  of  positive  and  negative  charges  depends 
on  the  dielectric  constant  and  the  pH  of  the  solvent.  The  interplay  of  long  range 
electrostatic  interaction  and  short  range  forces  of  other  origin  may  produce  sharp 
structural  transitions  upon  changing  certain  properties  of  the  solvent  [111-121]. 
The  interest  in  the  properties  of  PA  have  been  caused  not  only  by  their  rich 


structural  properties,  but  also  by  the  obvious  relevance  to  the  protein  folding 

problem  (122-128], 

The  overall  shape  of  a randomly  charged  linear  polyampholyte  chain  depends 
on  the  magnitude  of  the  net  charge  Q of  the  molecule  (6, 7],  If  this  charge  is  less 
than  the  critical  one  Qn  the  molecule  has  the  shape  of  a compact  globule,  where 
mutually  attractive  positively  and  negatively  charged  groups  nearly  compensate 
each  other.  The  globule  has  spherical  shape  in  order  to  minimize  its  external 
surface.  When  the  net  charge  exceeds  the  critical  value  Qr  ~ ZN,  the  globule 
splits  into  two.  This  new  configuration  has  larger  surface  area,  but  the  separa- 
tion of  two  resulting  globules  reduces  the  energy  of  electrostatic  repulsion.  This 
splitting  is  analogous  to  the  classical  instability  of  a charged  droplet  described 
by  1:  ieigh  (108].  Because  both  smaller  globules  arc  made  of  the  single  polymer 
chain,  they  can  not  be  too  far  apart  from  each  other  and  remain  connected  with 
a string.  At  larger  charges  each  globule  can  again  split,  therefore  forming  a linear 
chain  of  globules,  connected  strings.  The  number  of  globulae  is  controlled  by  the 
net  charge  Q.  The  optimal  number  of  globules  is  epproximately  [6,7] 

"= feVOa  (5.i) 

where  Qr  is  Rayleigh  charge  ( Qr  = 16 ir-yU3,  7:  surface  tension,  R:  radius  of  the 
globule ). 

The  randomness  in  the  distribution  of  charges  along  the  chain  results  in  fluc- 
tuations of  the  charge  in  each  globule.  These  fluctuations  (~  Nt12,  where  Nb 
is  the  number  of  monomers  in  a globule)  are  of  the  order  of  Qr  for  a globule. 
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Therefore,  the  fluctuations  caused  by  the  randomness  are  always  important  in  the 
size  distribution  of  globules-  There  is  another  factor  which  affects  the  statistics  of 
the  size  distribution  of  globules.  For  the  globules  on  a linear  chain,  the  effects  of 
electrostatic  repulsion  are  strongest  in  the  middle  of  the  chain,  and  weakest  near 
the  chain  ends.  As  a result,  the  size  of  globules  and  the  distance  between  them 
depend  on  the  position  of  these  globules  on  the  chain.  The  combined  effects  of 
randomness  and  position  dependence  make  any  analytical  and  numerical  study 
of  the  globule  distribution  in  the  linear  configurations  very  difficult  [6,7). 

5.2  Ring  Polymer 

In  this  work  I considered  the  configurational  properties  of  circular  chains. 
There  are  two  motivations.  First,  a ring  has  an  additional  symmetry  compared 
to  the  linear  chain,  since  all  monomers  have  the  same  structural  connections  to 
the  rest  of  the  chain.  The  effects  of  different  string  tension  discussed  above  can  be 
eliminated.  This  allows  the  study  of  the  asymptotic  size  distribution  of  globules 
on  smaller  chains  than  in  the  linear  case.  Second,  the  process  of  clusterization 
of  the  circular  polymer  might  provide  an  insight  into  the  micro-phase  separation 
in  polyampholyte  networks  (gels).  In  first  approximation,  these  networks  can  be 
thought  of  a collection  of  cross-linked  polyampholyte  rings.  In  the  same  way  as 
different  parts  of  the  ring  can  be  mutually  attractive,  different  chains  inside  the 
elementary  cell  of  a polymer  gel  network  can  be  attractive  and  form  a secondary 
network  of  compact  globules  connected  with  single  or  double  strings. 


First,  I will  neglect  effects  caused  by  disorder  and  assume  that  the  distribution 
of  charges  along  the  circular  chain  is  almost  homogeneous.  The  total  charge  of 
any  fragment  of  a polymer  of  length  m is  mQ/N,  This  is  similar  to  the  case  of  a 
homogeneously  charged  polyelectrolyte  in  a bad  solvent  (5),  where  the  monomers 
are  attracted  to  each  other  due  to  short  range  interactions,  but  repelled  due  to 
long  range  electrostatic  interactions.  The  only  difference  to  [5]  is  the  circular  form 
of  the  chain.  The  configuration  of  a polymer  is  controlled  by  the  overall  charge. 
If  this  charge  is  small,  the  polymer  is  in  a compact  globular  configuration,  if  the 
charge  is  large,  it  should  be  in  an  extended  multi-globule  configuration.  There 
are  two  major  possibilities:  One  possibility  is  that  the  extended  configuration  is 
linear,  similar  to  the  configurations  of  a linear  chain  (sec  Fig.  5.1  a — > V -*  c*). 
Another  possibility  is  that  the  multi-globule  configurations  will  have  a circular 
necklace  shape  as  shown  on  the  Fig.  5.1  a — t b — * c.  The  difference  in  the  energy 
of  these  configurations  depends  on  the  microscopic  details  of  the  model. 

Consider  the  energy  of  the  three  globule  configurations  of  a linear  and  a 
circular  polymer  (see  Fig.  5.1  c and  o').  There  are  three  major  contributions 
to  the  energy  of  these  configurations:  The  surface  energy  of  the  globules,  the 
surface  energy  of  strings  and  the  Coulomb  energy  of  the  interactions  between  the 
globules.  The  surface  energies  of  globules  of  both  configurations  arc  about  the 
same  and  I compare  only  the  two  remaining  contributions: 
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Figure  54:  Illustration  of  the  transitions  of  a ring  polymer.  When  the  charge 
of  the  polymer  increases,  new  configurations  can  be  formed  following  linear  of 
circular  sequences.  If  randomness  is  present,  there  is  a mixture  of  linear  and 
circular  patterns. 


Here  Q/ 3 is  the  charge  of  the  globules,  R is  the  distance  between  them,  (7;  (at) 
is  the  linear  tension  of  a double  folded  (single)  string.  This  tension  is  the  energy 
needed  to  pull  a unit  segment  of  the  string  out  of  the  globule.  Minimizing  both 
energies  with  respect  to  R,  I find 

El,„,„=2-/lQy/ai,  Bctrailar  = 0Qy/O\<  (5.3) 

If  the  double  folded  string  tension  at  satisfies  the  following  inequality, 

at  > 9/5 a,  (5.4) 

the  circular  configuration  has  lower  energy  compared  to  the  linear  configuration. 
If  the  reverse  inequality  holds,  the  linear  configuration  is  dominant.  The  ratio 
9/5  is  special  to  the  three  globule  configuration.  If  the  number  of  globule  n is 
large  (n  ^ 1),  the  energies  of  linear  and  circular  configurations  are 

Elmar  = » Log(tl)  + (n  - 1)0, fi,  Ea,alar  = n ioj(n)  + JUJ, . 

(5.5) 

In  this  case,  I have  the  inequality  at  > f7,  as  the  condition  of  lower  energy 
of  the  circular  configuration.  As  n increases,  the  energy  difference  6E  between 
linear  and  circular  configurations  converges  to  c(n)Q,  where  c(n)  is  a function  of 


c(n)  = 1{-/ai  - /ai)(Log(n))' 


(5.6) 


I expect  that  the  maximum  energy  needed  to  pull  out  a double  folded  string 
is  about  twice  the  energy  needed  to  pull  a single  string:  oj  ~ 2<T|,  the  actual 
ratio  depends  on  the  details  of  the  short  range  interactions  of  monomers  on  the 
two  strings.  If  the  monomers  are  interacting  through  saturatiug  forces  {with 
formation  of  H-bonds),  it  might  happen  that  the  unit  energy  of  the  double  string 
02  can  be  even  less  than  that  of  a single  string  0\- 

Summarizing,  for  charges  Q > Qr  the  minimum  energy  configuration  can 
be  linear  or  circular,  depending  on  the  details  of  the  model.  The  most  likely 
scenario  is  that  with  small  number  of  globule  n the  preferable  configuration  is 
linear,  while  for  large  n it  is  circular.  In  the  limit  case  when  the  electrostatic 
repulsion  between  monomers  dominates  over  all  other  interactions,  the  minimum 
energy  configuration  is  always  an  open  ring, 

5.4  Effect  of  a Random  Charge  Distribution,  Phase  Diagram 

Randomness  in  the  charge  distribution  creates  a fluctuation  of  the  net  charge 
in  each  globule.  Consider  a n globule  configuration.  The  typical  charge  fluctua- 
tion in  each  globule  is  ~ q,(Nln)'n.  The  fluctuation  in  electrostatic  energy  per 
globule  with  radius  r ~ (N/n)l^a  is  of  the  order  of  6Efiuct  ~ ql/&(N/n)2P.  If 
this  energy  SE/iMt  is  comparable  to  the  difference  between  the  energy  per  globule 
in  linear  or  circular  configuration  6E/n  = { — Eamiar)/n,  the  effects  of 
randomness  will  completely  control  the  local  arrangement  of  globules.  Setting 


6E/ i„e,  equal  to  SE/n,  I have 


c(n)0/n»»2/a(iV/n)1'3  (5.7) 

where  n ~ Q2/.\:.  This  equation  represents  a line  in  the  phase  diagram  plotted 
in  coordinates  Q2/Q7R  (=  n)  and  JV.  In  the  n 1 limit, 

IV*c(n)V/<&  (5.6) 

with  c(n)  given  by  Eq.  5.6.  The  Eq.  5.7  provides  an  upper  bound  for  the  region 
in  the  phase  diagram  where  the  effects  of  randomness  are  dominant  (Fig.  5.2). 
Above  this  line  (region  1)  the  minimum  energy  configuration  is  circular  or  linear, 
depending  on  the  microscopic  details  of  the  model.  The  dashed  vertical  lines  in 
the  region  I of  the  phase  digram  separate  regions  with  different  optimal  numbers 
of  globules  n.  In  this  region  there  are  also  numerous  meta-stable  configuration  of 
various  shapes.  The  energy  difference  between  the  different  topological  configu- 
rations is  of  the  order  ~ SE/n.  Below  the  line  given  by  Eq.  5.7  (region  II),  the 
ground  state  energy  configurations  are  completely  random.  In  addition,  there  are 
many  meta-stable  states  with  similar  energies.  The  lower  bound  of  region  II  is 
given  by  the  equation 

C?2  = c,N.  (5.9) 

Here  C2  is  a model  dependent  constant.  Below  this  line  (region  III)  the  charge  on 
the  chain  is  so  large  that  the  only  stable  configurations  are  open  rings  (see  Fig. 
5.2). 
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Figure  5.2:  Phase  diagram  of  a circular  polymer.  In  regions  I and  II  meta-stable 
states  are  observed.  In  region  I,  the  ground  state  is  either  linear  or  circular,  de- 
pending on  the  model.  In  region  II,  the  ground  state  is  controlled  by  randomness. 
In  region  III  (strong  net  charge  limit  Q > Ci N),  only  open  ring  configurations 
are  observed. 


5.5  Metastable  Configurations  and  Energy  Ba 


In  regions  I and  II  of  the  phase  diagram  there  are  many  configurations  which 
do  not  have  the  lowest  energy  but  are  relatively  stable.  Some  of  them  can  be 
relatively  easy  transformed  into  another  (soft  mode).  The  simplest  way  to  obtain 
a slightly  different  configuration  is  to  make  cyclic  shift  of  the  chain  between  each 
pair  of  nearest  globules  by  one  or  two  monomers.  This  transformation  preserves 
the  topology  of  the  overall  structure  of  globules  connected  by  strings.  The  differ- 
ence in  energy  of  the  new  configuration  is  of  the  order  of  ~ <to(Q/n)/«.  There  are 
about  N such  states  for  each  topologically  distinct  configuration.  Another  type 
of  transformation  is  related  to  a change  in  topology.  The  typical  barrier  between 
two  such  configurations  can  be  estimated  from  the  energy  needed  to  bring  two 
globules  together  and  then  regrouping  them  (see  Fig.  5.1).  This  energy  is  of 
order  £*  ■—  ( Q/n)7/a{N/n )*^3.  An  estimate  of  the  number  of  configurations  with 
different  topology  can  be  made  by  evaluating  the  number  N of  connected  clusters 
made  from  n globules.  This  is  the  well  known  “animals”  problem  [129],  and  the 
number  of  configurations  is:  N — n~°eAn,  where  the  exponent  a = 3/2  for  three 
dimensions  and  A is  a constant  of  the  order  of  one, 

5.6  Hysteresis 

If  the  overall  charge  of  the  chain  is  increased  (due  to  changes  in  solvent  proper- 
ties), the  formation  of  linear  structures  is  more  preferable.  Consider  the  transition 
from  two  to  three  globules.  This  transition  might  occur  by  splitting  one  of  the 
globules  with  formation  of  circular  or  linear  three  globule  configuration  as  shown 
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in  Fig.  5.1  c and  d.  Because  of  the  electrostatic  repulsion  forces  between  the 
globules,  the  splitting  of  the  globule  iu  parallel  direction  is  dynamically  more 
preferable  (Fig.  5.1  a -+  V ) than  the  splitting  in  perpendicular  direction  (Fig. 
5.1  a -+  5). 

The  linear  configurations  formed  in  this  way  may  not  be  the  ground  states. 
However,  they  are  meta-stable  and  dynamically  preferable.  The  time  required  to 
overcome  the  energy  barrier  between  linear  and  circular  configurations  is  propor- 
tional to  eEtlK‘T.  Depending  on  the  temperature  and  the  size  of  the  globules,  the 
transition  from  one  topology  to  another  cannot  be  observed  during  reasonable 

In  the  reverse  situation,  when  the  total  charge  decreases  from  a large  value, 
the  formation  of  a circular  configuration  is  dynamically  more  preferable. 

5.7  Monte  Carlo  Simulation 

My  model  and  the  details  of  the  Monte  Carlo  simulation  are  very  similar  to 
that  of  Kantor  and  Kardar  [6,7]  for  a linear  chain.  I consider  a ring  polymer  as 
a flexible  chain  of  IV  randomly  charged  monomers.  Each  monomer  has  charge 
+qo  with  probability  p and  charge  —go  with  probability  1 — p independent  of 
all  the  others.  A typical  random  sequence  of  N monomers  has  a total  excessive 
charge  of  about  (2 p - l)qo.  I characterize  the  sample  system  by  chain  length 
A' , temperature  T and  net  charge  Q.  At  a given  net  charge,  several  quenched 
configurations  are  tested.  The  interaction  between  a monomer  pair  i,j  is  the 
electrostatic  interaction  Ui-j  = 4,7,/r,  J During  the  simulation  the  monomers 


placed  on  a cubic  lattice.  Each  site  can  be  occupied  by  only  . 


and  the  length  of  the  link  between  two  subsequent  monomers  was  allowed  to  vary 
from  1 to  4 lattice  units.  The  Monte  Carlo  simulation  creates  new  configuration 
by  attempting  to  move  each  monomer  in  a randomly  created  sequence.  The  new 
configuration  is  accepted  with  probability  e~^E/koT,  where  A E is  the  energy 
difference  between  the  old  and  new  configurations. 

1 run  the  simulation  in  the  temperature  range  0.1  --  1 measured  in  q*/akg 
units  and  chain  sizes  .V  from  16  to  128.  Most  of  the  simulations  were  done  at 
temperature  T = 0.1  with  N = 96.  The  temperature  is  chosen  in  such  a way 
that  the  system  arrives  at  the  equilibrium  state  in  a reasonable  MC  time,  but 
low  enough  to  see  clear  transition  patterns. 

The specific  transition  patterns  are  different  for  the  various  models.  For  exam- 
ple, the  more  I increase  the  maximum  allowed  distance  between  the  subsequent 
monomer,  the  more  linear  configuration  will  be  preferable. 

I found  a very  rich  configurational  diversity.  This  is  caused  by  the  random- 
ness in  the  charge  distribution.  In  a lineal  polymer,  the  randomness  in  the  charge 
distribution  contributes  only  to  the  size  distribution  of  the  globules.  In  a ring 
polymer,  there  are  additional  possibilities  to  combine  two  remote  parts  of  the 
chain  containing  complementary  charged  monomers.  The  transitions  have  dif- 
ferent patterns  for  various  polymer  sizes.  As  a rule,  the  effect  of  randomness  in 
configurations  and  transitions  with  the  same  number  of  globules  n is  more  pro- 
nounced for  small  N.  For  N < 102,  the  transition  of  splitting  one  globule  was 
highly  irregular  for  different  sequences,  but  for  N > 102  the  transition  occurs  at 
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Figure  5.3:  Different  meta-stable  configurations  out  of  the  same  chain  with  fixed 
random  charge  distribution. 

Q « Qr  = VN,  as  predicted  by  mean  field  theory.  For  the  maximum  length 
N = 128,  the  2-globulc  to  3-globule  transition  was  affected  by  randomness  in 
the  charge  distribution  (see  phase  diagram  Fig.  5.2).  If  I increase  the  excess 
charge  at  fixed  Af,  the  number  of  globules  increases.  In  this  case  each  globule 
contains  less  monomers,  and  I found  very  irregular  cluster  shapes.  Some  sample 
configurations  are  shown  in  Fig.  5.3.  For  each  particular  charge  distribution,  dif- 
ferent meta-stable  configurations  can  be  accessed  by  repetition  of  the  cooling  and 
heating  procedure  described  below.  The  easiest  way  to  distinguish  different  con- 
figurations is  by  direct  observation.  It  is  very  difficult  to  assign  any  quantitative 
parameters  for  the  different  configurations. 
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Figure  5.4:  Plot  of  R*  vs.  Monte  Carlo  time.  Iu  the  beginning  of  each  cycle  for 
the  first  10s  steps  the  temperature  is  increased  to  T = 10.  Then  the  system  is 
at  T = 0.1  for  the  rest  of  the  cycle.  has  plateaus  at  various  levels  indicating 
different  meta-stable  configurations  shown  below  the  graph. 
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In  the  Fig.  5.4,  I plot  the  4th  moment  of  the  polymer  size  < flj  >s 
jfj  - Xj)*  while  the  system  undergoes  the  heating  and  cooling  cycles. 

In  the  beginning  of  each  interval  (1/5  of  the  interval),  the  system  is  heated  to 
T=10  and  configurations  become  nearly  Gaussian.  (<  Rj  > goes  down  for  these 
configurations!)  Then  I drop  the  temperature  to  T = 0.1  to  ensure  the  relaxation 
to  the  local  energy  minimum.  One  can  distinguish  the  configurations  by  the  value 
of  fij.  There  are  large  scale  plateaus  which  include  small  scale  fluctuations. 

The  linear  stretched  configuration  has  the  largest  Rj  (shown  as  the  highest 
plateau),  the  circular  open  shape  has  smaller  (middle  plateau).  The  config- 
uration with  a large  collapsed  part  or  a Gaussian  configuration  has  the  smallest 
value  of  R'9.  After  each  heating  and  cooling  process,  I reach  a different  configura- 
tion and  a different  level  of  plateau.  The  large  scale  jumps  imply  the  change  from 
a circular  to  a linear  configuration  or  vice  versa.  The  small  fluctuations  are  due 
to  the  shift  in  monomer  level.  With  chain  length  N = 64,  the  longest  relaxation 
time  f observe  to  reach  to  the  local  minimum  is  of  order  :V7  as  3000  — 4000  time 

Upon  increasing  the  net  charge  (Q  > Qj),  1 observe  a ring  configuration  with 
different  sizes  of  globules  as  a dominant  configuration.  I define  Q;  as  the  minimum 
net  charge  which  allows  exclusively  ring  configuration.  Qi  is  proportional  to  the 
polymer  size  N,  Qt  = cjsV.  The  model  dependent  constant  c2  is  approximately 
0.6  in  my  model. 


I found  that  ring  polyampholytes  have  remarkable  rich  configurational  prop- 
erties  as  a result  of  competition  between  tendencies  to  form  linear  or  circular 
structures.  The  different  parts  of  a ring  may  be  attracted  to  each  other  forming 
a compact  globule,  so  that  the  configuration  has  the  form  of  a random  cluster  of 
such  globules.  There  are  many  different  meta-stable  configurations  of  this  type 
for  the  same  chain  and  their  energies  are  about  equal.  I find  the  region  of  the 
phase  diagram  where  the  effects  of  randomness  in  the  distribution  of  charges 
along  the  chain  are  most  important.  I suggest  that  similar  multiple  conforma- 
tion phenomena  might  be  observed  in  a polyampholyte  gel.  In  the  same  way  as 
different  parts  of  a circular  chain  attract,  producing  various  configurations,  the 
different  parts  of  the  strands  inside  a gel  can  attract.  Local  rearrangement  of  such 
configurations  may  involve  strands  from  adjacent  elementary  cells  and  may  be 
thermodynamically  unlikely.  I expect  that  under  the  change  of  the  external  condi- 
tions (temperature,  pH,  salt  concentration,  etc.)  the  rearrangement  will  happen 
in  abrupt  steps  which  could  involve  the  whole  sample.  Nevertheless,  my  at- 
tempt to  link  the  multi-level  volume  phase  transition  observed  in  polyampholyte 
gels  [121]  using  the  analogy  with  single  chain  phenomena  should  be  considered 
with  caution.  In  my  approach  I neglected  the  possibility  of  screening  effects  due 
to  the  presence  of  counterions,  ions  of  salt  and  other  polyampholyte  molecules. 
Future  development  should  incorporate  these  effects  of  screened  interaction  as 
well  as  the  thermodynamical  contributions  of  other  ions. 


CHAPTER  6 

CONDENSATION  THEORY  ON  COUNTERIONS 


In  this  chapter,  I study  the  condensation  theory  of  counterions.  In  a polar  sol- 
vent, poiyelectrolytes  dissociate  into  charged  molecules  and  counterions.  When 
the  polyelectrolyte  chain  is  strongly  charged,  the  electrostatic  potential  near  the 
chain  is  large  and  some  of  the  counterions  remain  bound.  This  phenomenon  is 
known  as  “counterion  condensation  or  "Manning  condensation”  [8].  The  con- 
densation is  important  even  in  a very  dilute  solution  because  in  the  very  vicinity 
of  the  chain,  the  electrostatic  energy  of  counterions  (measured  in  kuT  unit)  can 
always  be  larger  than  unity.  Treating  conterion  condensation  analytically  is  dif- 
ficult because  it  occurs  in  the  regime  of  strong  Coulomb  interaction  where  the 
Debye-Hiickel  (linear)  approximation  is  not  valid. 

First  I review  the  well  known  solution  of  the  nonlinear  Poisson-Boltzmann 
(PB)  equation  in  a 2-D  cylindrical  cell  geometry  [130].  Then  I present  a new 
solution  of  the  PB  equation  with  different  boundary  conditions  appropriate  for 
the  isolated  polyelectrolyte  chain.  This  will  provide  a simple  explanation  of  a 
critical  charge  which  is  3/2  times  larger  than  in  Manning's  theory.  Finally,  I will 
discuss  how  the  interactions  between  two  test  charges  in  the  nonlinear  regime  are 
modified  due  to  the  presence  of  counterions.  This  result  is  then  related  to  the 
local  structure  (stiffness)  of  the  polyelectrolyte  chain. 


I consider  the  electric  field  produced  by  a long  rod-like  charged  polymer  in 
the  presence  of  counterions.  The  rod  consists  of  N monomers  with  radius  a,  and 
every  .4th  monomer  carries  unit  charge  e.  The  charge  parameter  A is 


A = Ig/aA 


(6.1) 


where  lg  = e*/[€kgT)  is  the  Bjerrum  length  which  is  a parameter  describing  the 
strength  of  electrostatic  interactions.  Bor  typical  polyelectrolytes,  the  Bjerrum 
length  is  about  7 A in  water. 

If  I ignore  .the  contribution  of  counterions,  the  electrostatic  potential  satis- 
fies Laplace  equation.  The  solution  is  obtained  by  Gauss  theorem  as  $(p)  = 
— 2A  log (,;//?,,).  The  counterion  density  around  the  chain  is  given  by  Boltzmann 
statistics  as  c(p)  = c0  exp(— $)  — p“\  The  total  number  of  counter  ions  per  unit 
length  within  a distance  p from  the  chain  is 


(6.2) 


When  the  line  charge  density  A is  larger  than  one,  the  above  integral  diverges 
in  the  vicinity  of  polymer  indicating  the  strong  condensation  of  counterions  on 
the  charged  rod.  As  the  condensation  proceeds,  the  counterion  clouds  effectively 
reduces  the  strength  of  Coulomb  interactions.  The  effective  value  of  the  charge 
seen  far  from  the  rod  is  loner  than  the  nominal  charge  density  on  the  polymer. 
The  condensation  stops  when  the  effective  line  charge  density  is  equal  to  unity. 


This  simple  condensation  theory  by  Manning  [8]  describes  the  condensation  as  a 
transition  between  two  states,  the  bound  state  and  the  free  state,  which  resets 
the  effective  charge  to  unity  for  large  charge  density.  Manning's  theory  shows  fair 
agreements  wnli  many  experiment.  However,  it  is  oversimplified  picture  [132, 133). 

in  a more  accurate  approach,  the  electrostatic  potential  of  the  polyelectrolyte 
chain  in  the  presence  of  counterions  should  be  described  by  the  PB  equation  (The 
potential  is  measured  iu  unit  of  kgT/e ). 

Vs$(p)  = -Mg  Y.  CilP)  (6.3) 

where  Cj  = c+(c_)  is  the  concentration  of  the  positive  (negative)  ions.  The  concen- 
tration is  determined  by  Boltzmann  statistics  c*(p)  = Cjcxp(— $(p)).  Therefore, 

V3$  = 4!re,//.-fl7'(c_  exp($)  - c+  exp(-4>))  (6.4) 

If  the  potential  is  small,  the  PB  equation  can  be  linearized,  resulting  in  the 
Debye-Hiickel  equation.  The  effective  potential  produced  by  a single  charge  is 
then  given  by 

*DB  = kBT^cxp(-Kp)  (6.5) 

The  contributions  from  different  charges  can  be  added  algebraically.  The  screen- 
ing length  k-1  in  related  to  the  Bjerrum  length  and  concentration  by  the  following 
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However,  in  the  vicinity  or  the  charged  polymer,  one  can  always  find  a region 
where  the  electrostatic  potential  is  large,  * > 1,  even  for  a very  dilute  solution. 
In  this  region,  the  nonlinear  aspect  of  PB  equation  becomes  important.  Thus, 
the  correct  condensation  theory  should  include  the  aspect  of  non  linearity.  In 
some  special  cases,  it  is  possible  to  solve  the  nonlinear  PB  equation  exactly  (130). 

6.2  Review  of  the  2-D  Cell  Model 

Fuoss  et  al.  have  studied  the  interaction  between  the  polyelectrolvte  chains 
and  counterions  using  the  PB  equation  [130,131],  They  assume  that  the  polymers 
are  infinite  rods  with  cylinder  radius  a and  are  all  parallel  to  each  other.  In 
between  these  rods,  there  must  be  equipotcntial  surfaces  where  the  electric  field 
vanishes(p  — /?,,}. 

Each  chain  with  its  surrounding  counterions  is  independent  of  all  the  other 
chains.  Therefore  the  problem  is  to  solve  the  electrostatic  problem  of  single 
rigid  cylinder  with  counterions  in  a 2-D  cylindrical  cell  [130].  The  electrostatic 
potential  in  this  cell  satisfies  the  following  Poisson-Boltzmann  equation. 


- -<cJexp(-$) 


with  the  boundary  condition 
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Figure  6.1:  The  geometry  of  the  cell  model.  The  electric  field  vanishes  at  the 
boundary.  The  polymer  is  approximated  as  an  infinite  rod  with  cylinder  radius 
a.  The  center  of  each  rod  is  located  at  the  center  of  each  cell.  The  net  charge 
wit*nn  each  cell  vanishes  due  to  the  electronical  neutrality  condition. 

Note  that  Fuoss  ct  al.’s  solution  is  controlled  by  two  parameters:  the  charge 
density  on  the  rod  A and  the  cell  size  Rt. 

The  two  solutions  below  and  above  Manning  condensation  are,  respectively, 

•t  = l°g[^^sinh2(tanh",^  + /3log(^))],  A<1  (6.10) 

*«  = lo8[^T  sinJ(tan-‘  P + /?log(|-))],  A > 1 (6.11) 

where  3 is  related  to  the  line  charge  density  parameter  A by  the  following  relation: 


The  distribution  of  counterions  can  be  calculated  using  Boltzmann  statistics  from 
the  potential  (6.10)  and  (6.11). 

These  solutions  indicate  that  there  is  no  bound  state  in  the  distribution  of 
counrerioiis.  However,  one  can  divide  space  into  two  regions:  1)  an  inner  region  in 
the  vicinity  of  the  chain  where  the  energy  is  larger  than  kgT,  therefore  counterions 
are  considered  as  bound,  and  2)  an  outer  region  where  the  interaction  is  smaller 
than  kgT,  therefore  the  counterions  are  moving  freely.  The  surface  potential 
$(p  = a),  (a  is  the  radius  of  the  cylindrical  rod)  increases  almost  linearly  with 
the  increase  of  the  line  charge  density  A until  A reaches  unity.  For  large  line 
charge  density  A > 1,  $(/>  = a)  increases  logarithmically. 

It  was  shown  that  when  0 < A < 1 the  effective  charge  density  increases  to 
unity  as  the  condensation  proceeds.  On  the  other  hand,  if  A > i , the  effective 
charge  density  is  limited  by  1/A  (46).  In  the  infinite  cell  model  (in  the  limit 
Rb  -*  oo)  [134],  Russel  also  showed  that  the  charge  density  is  limited  by  unity 
where  the  spacing  between  the  charges  is  the  Bjcrrum  length. 

6.3  New  Solution  of  the  Poisson-Boltzmann  Equation 

In  this  section,  I present  a new  solution  which  satisfies  different  boundary  con- 
ditions. I will  show  that  this  new  solution  describes  a different  physical  situation. 
I consider  a dilute  solution  of  polyelectrolytes  where  polyelectrolyte  chains  are 
isolated  from  each  other.  In  this  limit,  one  can  still  use  cylindrical  coordinates 
near  the  charged  rod.  But  in  a distance  comparable  to  the  chain  length,  the  3-0 
boundary  condition  must  be  considered.  The  PB  solution  in  the  vicinity  of  the 


Figure  6.2:  The  polyelectrolyte  chains  in  a dilute  solution.  The  electric  field  does 
not  vanish  at  the  boundary. 

cylinder  does  not  need  to  satisfy  the  charge  neutrality  condition  at  the  spherical 
boundary  shown  on  Fig.  6.2. 

I use  cylindrical  coordinates  p,  Q,  z in  the  vicinity  of  the  charged  rod-like 
polymer  (all  lengths  are  measured  in  terms  of  the  cylinder  radius  n).  The  po- 
tential outside  of  the  cylinder  is  described  by  the  PB  equation  (Eq.  6.4).  In  the 
vicinity  of  the  cylinder  (p  < pe  <£  d),  the  potential  is  large.  3>  1,  and  only  the 
first  term  of  the  r.h.s.  of  Eq.  6.4  is  important;  therefore, 

V2$,  = 4uec-  cxp(4>,)  (6.14) 

where  index  "i"  indicates  that  the  solution  is  valid  inside  of  the  cylindrical  regime 
(p  < pc  < d).  In  the  vicinity  of  the  charged  rod.  the  effect  of  finite  length  of  the 
polymer  is  negligible.  I found  the  solution  of  this  equation  for  the  infinite  length 


rod  has  the  very  simple  form. 


•W-Wl 5^  <«•“> 

The  solution  has  two  tuning  parameters  b and  C.  There  is  a simple  relation 
between  b and  the  charge  density  on  the  surface. 

*— mo 

The  parameter  b gives  the  location  of  the  hypothetical  surface  inside  of  the 
charged  cylinder  where  the  potential  diverges.  The  asymptotic  behavior  of  the 
electrostatic  potential  (p  ->  oo)  is  0 = 31og(l/p).  The  effective  charge  density 
seen  from  far  distance  of  the  cylinder  is  universal  value  AM  = 3/2.  The  counterion 
density  in  the  solvent  (p  > 1)  is  given  by 


(6.17) 


The  total  number  of  counterions  at  the  distance  p from  the  rod  can  be  calculated 
by  integrating  the  above  equation  Eq.  6.17:  c,(p)  = ff  c-[p)2npdp.  The  main 
contribution  to  this  integral  comes  from  the  region  near  the  charged  rod.  The 
solution  Eq.  6.15  is  valid  only  for  strongly  charged  polymers  with  charge  A > 3/2. 
When  the  charge  density  A < 3/2,  the  PB  equation  can  be  linearized  and  the 
solution  has  the  form  —A  log(p).  The  contribution  of  counterions  is  negligible  in 
this  dilute  solution  limit. 


6.4  3-D  Discussion 


The  constant  C in  9,  will  be  determined  by  the  second  boundary  condition 
which  relates  the  inside  solution  * to  the  solution  outside  of  the  boundary.  I 
consider  another  hypothetical  cylinder  characterised  by  the  radius  p,  and  length  d 
surrounding  the  charged  rod.  For  the  outside  of  this  cylinder(p  > pc,  |z|  > d/2), 
the  potential  is  expected  to  be  small  (4  < 1).  In  this  region,  the  counterion 
density  (Eq.  6.17)  is  very  small  so  that  I can  consider  the  Laplace  equation  of  a 
charged  rod  with  constant  line  charge  density  3/2.  The  electrostatic  potential  in 
this  regime  can  be  obtained  by  calculating  that  of  a charged  cylindrical  rod  with 
radius  pc  and  length  d.  Therefore,  the  electrostatic  potential  is 

= | (6-16) 

where  !•<(*>)  is  the  smalier(larger)  of  d/2  and  r. 

I match  the  boundary  condition  in  the  3-D  geometry  at  a finite  distance  from 
the  center  of  the  rod.  Therefore,  my  new  solution  describes  the  regime  where 
there  is  no  overlap  between  the  chains.  The  picture  of  counterion  condensation 
in  3-D  is  different  from  the  2-D  cell  model.  Since  the  counterion  has  to  be 
brought  from  3-D  space  rather  than  2-D,  the  entropy  loss  of  a counterion  in  3-D 
is  of  the  order  of  ~ 3 log(p)  instead  of  2 log(p)  in  2-D.  In  Manning  theory  the 
condensation  occurs  in  the  cylindrical  cell  geometry  when  the  electrostatic  energy 
gain  of  bringing  the  counterion  to  the  surface  2A  log  p becomes  comparable  with 


the  2D  entropy  loss  of  the 


s.  My  result 
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dilute  solutions  of  polyelectrolytes.  In  this  case  the  counterion  should  be  brought 
from  the  3D  space  and  the  entropy  loss  is  3A  logp  which  is  3/2  times  larger  than 
in  Manning's  case.  Therefore,  the  onset  of  condensation  occuts  at  higher  charge 
density. 

If  the  end  effects  are  not  neglected,  the  constant  C has  z dependence.  The 
electrostatic  potential  at  the  edge  diverges  logarithmically: 

Up  = a,i)  = \ log(i^)  + (6.19) 

where  p = 


6.5  Interaction  between  Test  Charges 

I will  show  that  due  to  the  enhanced  concentration  of  ions  near  the  chain 
the  actual  interaction  is  different  from  Debye-Huckel  type  interaction.  Consider 
a small  test  charge  on  the  charged  cylinder.  Due  to  the  presence  of  the  test 
charge  there  is  a small  perturbation  in  the  potential  proportional  to  the  change 
of  charge  density.  The  perturbed  electrostatic  potential  can  be  written  as  a sum 
of  the  non-perturbed  potential  and  small  perturbation:*  = *„  + *,.  The  PB 
equation  for  $i  can  be  written  as: 

V2*,  = y exp(*„)*,  (6.20) 

The  interaction  of  charges  near  the  chain  are  controlled  by  two  new  charac- 
teristic scales.  One  of  them  is  the  size  of  the  nonlinear  region  pc  near  the  chain 
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where  the  Dcbye-Huckel  approximation  is  not  valid.  Another  scale  is  the  “local” 
Debye  screening  length  R‘d  near  the  chain.  By  comparison  with  Eq.  6.4,  I can 
define  the  local  Debye  length  by 

75  = (*J/2ex  p(*.))-''J  (6.21) 

where  is  the  non-perturbed  potential.  Therefore  the  interaction  between  two 
test  charges  must  be  modified  by  using  the  local  Debye  length.  If  the  distance 
r between  two  test  charges  is  less  than  Rj,  the  interaction  between  them  is  not 
screened. 

By  taking  the  derivative  of  Eq.  6.15, 1 obtain  the  change  of  the  potential  due 
to  the  small  perturbation  in  charge  density. 

<f<fi  ~ 5b/ p ~ 5X/p  (6.22) 

This  relation  shows  that  the  screening  of  counterions  reduces  the  interaction 
strength  by  factor  of  1/p.  Suppose  one  of  the  charges  is  on  the  chain  and  another 
is  placed  at  a distance  r ( r > flj)  away  in  perpendicular  direction.  The  field 
line  starting  at  the  charge  on  the  cylinder  can  reach  the  other  charge  with  a 
probability  1/r.  Due  to  the  counterion  condensation,  the  strength  of  the  field 
lines  is  reduced  by  an  additional  factor  1/r  given  by  Eq.  6.22.  Therefore,  the 
interaction  decays  as  1/r’.  If  the  second  charge  is  placed  on  the  cylinder,  the  field 
lines  passing  near  the  cylinder  have  little  chance  to  get  to  the  other  charge.  Two 
test  charges  on  the  cylinder  are  interacting  through  field  lines  detouring  through 


the  outer  region  (r  > flj).  The  interactions  will  be  reduced  by  an  additional 
factor  of  1/r,  leading  to  an  overall  decay  of  1/r3. 


I confirmed  this  argument  by  a numerical  path  integral  method.  In  a com- 
puter simulation,  the  random  walk  starting  from  a point  on  the  cylinder  in  the 
electrostatic  field  6.15  has  a probability  1/r3  to  get  to  the  other  point  at  a dis- 
tance r along  the  perpendicular  direction.  The  probability  is  1/r3  if  the  second 
point  is  at  a distance  r on  the  cylinder.  If  the  distance  between  the  charges  is 
larger  than  pe  the  interaction  is  described  by  the  Debye-Hiickel  formula. 


CHAPTER  7 

PERSISTENCE  LENGTH  OF  A POLYELECTROLYTE  CHAIN 


I present  a first  principles  calculation  for  the  persistence  length  of  a linear 
chain  in  both  the  strongly  and  weakly  charged  case.  In  the  presence  of  salt  ions, 
using  curved  cylindrical  coordinates  we  calculate  the  change  in  the  electrostatic 
potential  and  ion  concentration  in  the  vicinity  of  a bent  chain.  I find  the  per- 
sistence length  to  dependence  on  the  Debye  length  as  I,  ~ it"1  for  a strongly 
charged  polymer.  This  result  can  be  explained  by  enhanced  screening  effects  due 
to  the  high  concentration  of  ions  in  the  vicinity  of  the  chain,  compared  to  the 
concentration  in  the  bulk.  For  a weakly  charged  polymer,  we  recover  the  OSF 
result  l,  - *-s  [81,82]. 


7.1  Persistence  Length 

The  difficulty  in  treating  charged  polymers  lies  in  the  noil-trivial  influence 
of  the  long  range  Coulomb  interaction  on  the  local  structure  of  polymer.  The 
persistence  length  lp  is  the  parameter  describing  the  flexibility  of  the  polymer 
chain  [1,3].  The  persistence  length  is  defined  as  the  decay  length  of  the  correlation 
length  of  the  tangential  vector.  The  persistence  length  can  be  deduced  from  the 
measurement  of  various  physical  quantities  [89-93].  Many  theoretical  works  have 
attempted  to  calculate  the  persistence  length  of  a charged  polymer  [81-84, 86- 
88, 131],  Under  the  presence  of  electrostatic  interactions,  the  bending  brings  two 
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of  the 


(7.1) 


Fixman  (82) 


contribution  lOSF,  i,  = fo  + W,  with 

losr  = (7-2) 


has  large  bare  persistence  length.  However,  for  intrinsically  flexible  polyelec- 
trolytes, however  shows  strong  deviation  from  OSF  theory.  As  one  decreases  the 
ion  concentration  c,  the  persistence  length  increases  as  c~'.  Using  the  variational 
method,  Barrat  and  Joanny  [36]  argued  that  for  intrinsically  flexible  chains,  the 
fluctuations  in  chain  configuration  can  have  substantial  influence  on  the  rigidity  of 
the  chain,  therefore  the  persistence  length  remains  of  the  order  of  the  electrostatic 
screening  length.  This  result  is  consistent  with  the  picture  of  an  electrostatic  blob 
suggested  by  de  Gennes  et  al  (see  Fig.4.1  a).  Inside  of  the  electrostatic  blob  the 
thermal  energy  is  comparable  to  the  Coulomb  interaction  so  that  the  chain  re- 
mains undistorted  by  electrostatic  interaction.  However,  the  calculations  by  Li 
and  Witten  [88]  show  that  the  effect  of  conformational  fluctuations  of  a flexible 
ring  polymer  does  not  change  the  OSF  result. 

On  the  other  hand,  several  computer  simulations  indicate  that  the  pertur- 
bation by  the  counterions  is  much  stronger  than  indicated  by  the  DH  formula 
[94-97].  This  implies  stronger  screening  of  the  electrostatic  interactions  and 
shorter  persistence  length.  More  recent  molecular  dynamics  simulations  by  Micha 
and  Kremer  [98, 99|  show  that  the  persistence  length  has  sublincar  dependence 
on  the  screening  length:  Ip  ~ K~v,y  < 1. 

7.2  Debye  Hiickel  Approximation 

It  is  clear  that  there  arc  limitations  in  applying  DH  theory  to  a polymer 
system.  The  Debye-Huckcl  potential  is  obtained  as  a solution  of  a linearized 


; of  the  Poisson  - Boltzmann  equation: 


= k2$  (7.4) 

In  the  actual  problem,  in  the  vicinity  of  the  chain  the  potential  $ is  larger  than 
ksT/e  and  the  PB  equation  becomes  nonlinear.  Therefore  in  the  vicinity  of 
the  chain  the  concentration  of  ions  is  higher  than  in  the  bulk.  The  screening 
effects  should  be  enhanced.  Le  Bret  [83)  and  Fixman  [84]  have  calculated  the 
persistence  length  numerically  using  the  nonlinear  Poisson-Boltzmann  equation 
rather  than  the  linearized  Debye-Huckel  approximation.  The  numerical  value  of 
the  persistence  length  shows  strong  deviation  from  the  OSF  result  at  high  ionic 
concentration  [84). 

In  the  following  sections,  I discuss  the  influence  of  the  nonlinearity  of  the  PB 
equation  on  the  stiffness  of  the  chain.  I extended  the  analytic  approach  to  obtain 
analytic  expressions  for  the  cross  over  from  "k"3”  (rigid  chain)  to  V1”  (flexible 
chain)  behavior  of  the  persistence  length.  I perform  a straightforward  calculation 
of  the  distribution  of  the  potential  and  the  ion  concentration  in  the  vicinity  of  a 
slightly  bent  charged  polymer  chain  with  curvature  radius  H.  Then,  I calculate 
the  bending  energy  by  expanding  the  PB  equation  up  to  the  second  order  of 
1/R  in  curved  cylindrical  coordinates.  This  expansion  procedure  is  similar  to 
Fixman's  approach  as  described  in  [84].  It  appears  that  the  main  contribution  to 
the  rigidity  comes  from  the  region  of  radius  r < pc  near  the  chain. 


Numerical  Solution  in  Cylindrical  Coordinates 


In  this  section,  I study  the  nonlinear  PB  equation  numerically  in  cylindrical 
geometry  in  order  to  illustrate  the  nonlinear  effects  of  the  PB  equation,  i consider 
a single  infinite  charged  rod  with  radius  a with  line  charge  density  A.  With  Eq.6.4, 
I write 

VJ$  = ysinh(<h)  (7.5) 

with  the  boundary  conditions 


9(p)  !,==.  = 0.  (7.7) 


Following  the  ideas  of  the  previous  chapter,  I confirmed  that  the  solution  of 
the  nonlinear  PB  equation  has  following  properties,  1)  When  A < 1 the  solution 
is  analytic  everywhere.  2)  When  A > 1 the  solution  has  singular  surface  inside 
of  the  cylinder  at  the  distance  r = 6 from  the  cylinder  center.  The  parameter  b 
has  also  weak  dependence.  The  location  of  the  singular  surface  b is  related 
to  a-1  via  the  following  equation: 


if*-1)  - = oo)  a a*/*-1.  (7.8) 

There  is  a difference  to  the  cell  model  where  the  solution  is  controlled  by  two 
independent  parameter  A and  Rt,.  I point  out  that  this  solution  of  the  non-linear 
Poisson  Boltzmann  equation  is  controlled  by  only  one  parameter:  6/k~'.  3)  If 


Figure  7.1:  a)  The  solution  of  the  PB  equation  for  A = 1,1.2, 1.4, ..  increasing 
from  the  lower  to  the  upper  Fig.  b)  The  counterion  distribution  as  a function  of 
distance  from  the  cylinder  for  various  values  of  A.  The  concentration  decreases 
dramatically  in  the  vicinity  of  cylinder.  The  counterion  density  in  the  region 
away  from  the  cylinder  is  not  affected  by  the  charge  density  change. 


b/n~l  « 1,  the  solution  of  the  nonlinear  PB  equation  has  universal  properties. 
Far  away  from  the  rod,  the  solutions  with  different  A are  on  top  of  each  other.  As 
A increases  the  solutions  are  emerging  to  a universal  curve  (see  graph  7.1).  The 
two  graphs  7.1  a)  and  b)  show  the  counterion  density  and  potential  distribution 
for  various  value  of  surface  charge  density.  The  influence  of  large  A can  be  seen 
only  near  the  cylinder.  The  potential  near  the  cylinder  is  much  larger  than  A-pT 
and  counterions  density  is  strongly  enhanced  in  this  regime.  The  change  of  elec- 
trostatic potential  and  counterion  density  upon  the  change  of  the  surface  charge 
density  is  negligible  at  a distance.  As  predicted  by  other  theoretical  approaches, 
we  can  clearly  see  that  the  redistribution  of  counterions  compensates  the  excess 
charges. 


7.4  Calculation  of  the  Bending  Energy 

In  order  to  obtain  the  persistence  length,  I follow  Fixman's  expansion  method 
[84]  using  the  curved  cylindrical  coordinates  with  radius  of  curvature  R.  The 
geometry  of  curved  cylindrical  coordinates  is  shown  on  the  Fig.  7.2.  The  PB 
equation  in  this  curved  cylindrical  coordinates  can  be  written  in  the  following 

p(  1 — pacosip)dp^  pacostjt)  4-  p(l  — pa  cos  0)  d<t>p(  1 — /XI  cos  o)  do 
(7.9) 


= 4trep_exp('h_) 


Figure  7.2:  The  geometry  of  curved  cylindrical  coordinates  are  shown.  The 
curvature  a is  the  inverse  of  the  radius  R. 

The  electrostatic  potential  is  expanded  witli  respected  to  the  inverse  of  the 
radius  of  curvature,  a = 1/R. 

#(*>)  = Y.  *»"•(/’)  ros( (710) 

Here  the  electrostatic  potential  $ is  measured  in  terms  of  kgT/e.  All  lengths  will 
be  measured  in  terms  of  the  radius  of  cylinder  a.  Assuming  that  the  bending 
curvature  is  very  small  a<l,  we  expand  $ up  to  the  second  order  of  a and 
neglect  the  higher  order  of  a: 


Mp)  = *o(p)  + *11  (p)  cos  0a  + ($m(p)  + (/»)  cos  2d)a2 


(7.11) 


The  electrostatic  contribution  to  the  persistence  length  is  obtained  from  the  en- 
ergy change, 

l,  = ^A  E/k„T.  (T.12) 

The  total  electrostatic  energy  of  the  ion  system  in  a given  state  is  the  work 
required  to  raise  its  charge  density  from  zero  to  its  final  value  A.  The  electro- 
static energy  increase  per  unit  length  due  to  the  bending  can  be  obtained  by  the 
electrostatic  energy  difference  between  the  straight  (40  and  bent  ($')  cylinders. 

AE  = e|dn|VUi-  ®|p=i  )dA  (7.13) 


where  /dfl  is  the  integration  over  the  azimuthal  angle.  After  angle  integration, 
all  odd  functions  will  vanish.  Therefore,  the  electiostatic  contribution  to  the 
persistence  length  comes  from  the  second  order  term.  The  energy  of  bending  is 

A E/kgT  = ej  dfl^VaoU,  + ^?ile=i)dA  (7.14) 


To  first  order  of  a cos  ^ we  have 


OTo  1 a-t-n 
dp  + p dp 


(7.15) 


The  solution  of  Eq.7.15  is  almost  independent  of  A in  a distance  far  from  the 
rod  (approximate  universality).  To  second  order  in  a.  I consider  all  isotropic(m  = 
0)  terms. 


(7.16) 


where  the  source  term  S[p)  is  given  by 


First,  let  me  consider  n weakly  charged  polyelectrolyte  (A  < 1).  The  last  term  of 
Eq.7.16  is  negligible  because  there  is  no  counterion  condensation,  therefore  the 
screening  by  the  counterions  is  weak.  I assume  that  there  is  no  significant  charge 
redistribution  upon  bending,  this  leads  the  initial  condition  The 
source  term  S(p)  is  approximately  ~ A.  Therefore  Eq.  7.16  can  be  approximated 


tFfra,  1 _ 
d?  + ~p  dp  ~ 


(7.18) 


The  last  term  of  Eq.  7.16  Becomes  important  when  <I>M  as  1/y/H.  This  means 
that  the  amplitude  of  should  vanish  at  the  screening  length  p a *-1.  By 
integrating  Eq.  7.18,  we  obtain  the  surface  potential  $20 [p  — 1)  a Arc'2.  In  order 
to  obtain  the  persistence  length,  we  then  integrate  with  respect  to  A,  /0‘  <I>MrfA. 
I recover  the  OSF  result  for  weakly  charged  polymer 


(7.19) 


For  the  strongly  charged  polyelectrolyte  (A  > 1),  the  interaction  range  s'1 
is  replaced  by  the  local  Debye  radius  (re*)-1,  which  is  defined  by  Eq.  6.21.  The 
interaction  range  (s')'1  is  dramatically  reduced  by  the  increased  concentration 
of  counter  ions  near  the  cylinder  where  the  electrostatic  potential  is  large  (4>  » 
l.p  < (x*)-1).  1 can  see  that  there  are  two  physically  different  regions.  In  the 


first  regime  (the  very  vicinity  of  the  cylinder),  the  PB  equation  consists  of  two 
dominant  contributions 


+ — - y «**®2o  = 0,  «5,U, 


(7.20) 


The  second  term  («')a  = Ye>°^  >s  operant  only  in  the  very  vicinity  of  the 
cylinder  p < (sc*)-1.  As  a result  of  the  integration,  we  have 


(7.21) 


where  we  use  the  identity  yc^’pdp  = (A  - 1).  Outside  of  this  region,  the 
.situation  is  same  as  in  the  weakly  charged  case.  As  in  the  case  of  Manning 
condensation,  the  charge  density  A is  replaced  by  unity. 


Same  as  in  the  weekly  charged  poiymer,  $20  is  of  the  order  of  unity  at  p = 
Suppose  we  approach  the  charged  cylinder  from  the  distance.  The  electrostatic 
potential  is  increasing  monotonically  until  p aj  (A  - l)d>2o  As  we  get  closer  to 
the  cylinder,  the  screening  term  becomes  important  so  that  the  magnitude  of  the 
potential  starts  decreasing.  I can  replace  by  a constant  and  solve  the 

differential  equation.  In  the  weakly  charged  case,  the  potential  should  be  of  the 
order  of  unity  at  p = 1.  Therefore,  we  find  the  value  of  $20  on  the  cylindrical 
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Figure  7 3:  Persistence  length  vs.  Debye  length  as  given  by  Ec.  7.23  with  A = 4.2. 
The  data  points  are  taken  from  Fixman’s  numerical  calculation. 

In  order  to  obtain  the  persistence  length  for  the  strongly  charge  rod  (A  > 1), 
we  have  to  integrate  ®f0  and  $5o  for  the  0 < A'  < 1 and  1 < A'  < A regions 
respectively.  The  persistence  length  is  the  linear  combination  of  the  two  integrals, 
we  obtain  following  formula, 


Ip  = <«-»  + bK-'  log(A-l).  (7.23) 

I fit  the  parameters  a,  b by  the  comparison  with  Fixman’s  numerical  results  for 
two  different  A.  I obtain  a as  0.604  and  b =s  7.23  which  allow  the  best  fitting  for 
both  of  A.  The  result  shows  good  agreement  with  Fixman’s  numerical  calculation 
(see  Fig.  7.3).  By  analysing  the  formula  Eq.  7.23,  I find  that  the  crossover 
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occurs  when  the  screening  length  is  approximately  k_1  ~ 10a.  For  a small  salt 
concentration  k~1  > 10a,  the  second  term  of  the  Eq.  7.23  is  dominant.  At  a 
large  salt  concenumion  n - < 10a  and  the  persistence  length  shows  a linear 
dependence  on  «“*. 

Summarizing,  using  a qualitative  analysis  of  the  perturbation  expansion  in 
curved  cylindrical  coordinates,  we  obtained  a simple  phenomenological  formula 
which  describes  the  crossover  of  the  persistence  length  from  the  weakly  to  the 
strongly  charged  limit. 


CHAPTER  8 
CONCLUSIONS 

In  this  thesis  I studied  the  conformational  properties  of  polymers.  I first 
discussed  two  polymer  systems  in  which  the  configurational  properties  arc  inti- 
mately related  to  the  dynamic  properties.  In  the  problem  of  a polymer  chain 
penetrating  a membrane  I calculated  the  probabilities  of  different  configurations 
of  the  polymer  and  the  corresponding  mobilities.  1 found  the  optimal  configura- 
tions which  give  the  main  contribution  to  the  transport  through  the  membrane 
in  both  Rouse  and  Zimin  dynamics.  This  process  of  diffusion  through  the  mem- 
brane can  be  used  for  separation  of  polymeis  with  different  length.  I showed  that 
the  total  permeability  has  a power  law  dependence  on  the  molecular  weight. 

I also  studied  the  conformations  of  DNA  molecules  undergoing  gel  electrophore- 
sis. In  a strong  electric  field,  the  configuration  of  DNA  consists  of  many  growing 
and  shrinking  branches.  The  power  law  behavior  of  the  branch  size  distribu- 
tion confirms  the  hypothesis  of  a self  similar  tree-like  structure.  It  is  observed 
that  DNA  has  a periodic  evolution  pattern  and  stays  in  a simple  uV"-shaped 
configuration  during  a considerable  fraction  of  the  period.  Although  the  mobil- 
ity of  DNA  is  independent  of  molecular  weight  in  constant  field  electrophoresis, 
the  knowledge  of  the  configuration  at  a given  time  will  help  to  develop  more 


ids  of  electrophoretic  separation  of  DNA. 


In  the  second  part  of  the  thesis,  I studied  the  conformational  properties  of 
charged  polymers.  I considered  a single  charged  polymer  which  is  isolated  from 
other  polymers. 

I showed  that  a polyampholyte  chain  has  a glassy  phase  where  the  config- 
urations are  controlled  by  the  randomness  in  the  charge  distribution  along  the 
chain.  The  metastable  configurations  in  this  phase  are  a result  of  the  competition 
between  the  tendencies  to  form  linear  or  circular  structures.  I suggested  that  sim- 
ilar multiple  conformation  phenomena  could  exist  in  polyampholyte  gels,  which 
could  be  connected  to  the  volume  phase  transitions  observed  in  the  experiment 
by  Tanaka  [121].  By  solving  the  Poisson-Boltzmann  equation,  I investigated  the 
phenomenon  of  counterion  condensation  and  the  resulting  stiffness  of  a poly- 
mer. The  flexibility  of  a strongly  charged  polyelectrolyte  can  be  explained  by 
the  strong  screening  of  the  Coulomb  repulsion  due  to  the  enhanced  counterion 
concentration  near  the  polymer.  I also  described  a new  condensation  scenario 
for  the  dilute  solution  in  which  the  critical  charge  is  3/2  times  larger  than  in  the 
classical  Manning  theory. 

The  understanding  of  the  microscopic  picture  of  a single  polymer  is  the  first 
step  towards  a study  of  melts  or  polymer  solutions  in  semidilute  or  concentrated 
regime,  which  is  essential  for  many  practical  purposes.  Many  properties  of  poly- 
electrolyte solutions  are  controversial.  Much  more  work  is  needed  to  reach  a level 
of  understanding  equivalent  to  that  obtained  for  neutral  polymers.  In  a polymer 
solution,  interchain  interactions  exist  as  well  as  intrachain  interaction.  The  study 
of  many-chain  effects  including  entanglement  effects  will  be  the  subject  of  future 


study.  In  addition  to  the  analytical  approach,  the  importance  of  numerical  stud- 
ies is  obvious.  The  numerical  simulation  can  include  counterions  explicitly  which 
goes  beyond  the  Debye-Huckel  or  Poisson-Boltzmann  approximation. 


APPENDIX  A 
PROOF  OF  EQUATION  2.1 


In  the  absence  of  excluded  volume  interactions  the  relation  pm  = \p,r/  can  be 
verified  by  the  following  simple  gedanken  experiment:  Let  me  assume  that  the 
same  polymer  concentration  is  established  on  both  sides  of  the  membrane,  but 
the  polymers  on  one  side  are  slightly  different  from  polymers  on  the  other  side 
of  the  membrane  (e.g.  they  can  be  deuterated).  The  polymers  of  both  types 
start  to  diffuse  and  mix  each  other.  The  total  concentration  of  both  types  of 
polymers  will  be  the  same  as  in  equilibrium,  but  the  concentration  profile  of 
each  component  can  be  calculated  only  if  the  dynamics  of  the  polymer  morion  is 
known.  In  the  middle  of  the  membrane  the  concentrations  of  each  components 
should  be  equal  to  each  other  because  of  the  symmetry  of  the  problem.  The  sum 
of  this  two  concentrations,  i.e.  the  total  concentration,  should  be  the  same  as  in 
equilibrium.  This  means  that  the  concentration  of  each  component  is  equal  to 
one  half  of  the  equilibrium  concentration  as  stated  in  Eq.  2.1. 
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APPENDIX  B 

CALCULATION  OF  PARTITION  FUNCTION 


When  the  polymer  crosses  the  membrane,  it  must  penetrate  the  membrane  at 
least  once.  I calculate  the  partition  function  Z,(p,Wi)  ( Z2(/s,  rV2) ) of  a polymer 
of  length  Ni{Ni)  attached  to  the  membrane.  Then  I write  the  partition  function 
of  the  whole  polymer  Z(p,  N)  as  a product  of  two  partition  functions: 


Z(p,iV)  = 2,(p,iV1)Z2(/i,Af3)e-' 


(B.1) 


One  may  set  Ni  = iVj  = iV/2,  in  which  case  the  constraint  from  the  membrane 
is  maximal.  The  weight  e~'1  is  attributed  to  the  single  membrane  crossing. 

From  the  well  known  solution  of  the  first  passage  problem  [135],  I calculate 
the  Green  function  of  random  walk  of  length  n which  starts  at  the  surface  of  a 
membrane  and  does  not  cross  the  membrane  again.  One  also  needs  the  Green 
function  of  a walk  which  starts  and  ends  on  the  membrane  and  does  not  have 
any  other  contacts  with  membrane  in  between.  I call  these  Green  functions  Ga 
and  G4,  respectively  (see  Fig.  B.l  a and  b): 


(B.2) 


The  configurations  contributing  to  the  partition  function  Zi(/t,  ;Vt)  are  shown  in 


Fig.  B.l  c. 


» were  studied  for  the  adsorption  problem  [20,21  j. 
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S r0-! 


Figure  B.l:  a),  b)  Graphic  representations  of  the  functions  G0,  Gt.  c)  Diagram 
series  of  all  contributions  to 

The  statistical  weight  of  each  of  them  can  be  written  as  a multiple  product  of  Ga 

Zxb,  «i)=  £ e"*''/c.(n0)C,(n1)...Gs(n,J4(rV1  - £n,)*l0dn1...dnll 

(B.3) 

When  kt  = kt  = 0,  the  polymer  passes  the  membrane  only  once.  The  total 
number  of  crossings  is  k,  +i2  + 1.  In  terms  of  Zt (p, Ni), Wj),  < 1/A  >=< 
l/(*i  + ki  + 1)  > is  given  by 

< 1/1  >=  jT  d/ndin,  tf(/a,  - w)  W,)  Zj(W,  Ar,)|Wi=1V]  (B.4) 

Let  me  define  the  Laplace  transform  Z,(p,  s)  of  the  one-side  partition  function. 
The  evaluation  of  the  partition  function  Z|(p,  s)  can  be  made  by  means  of  a 
diagram  technique  (see  Fig.  B.l  c).  The  Laplace  transforms  of  G„(p,n),Gs(/i,n) 
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(B.5) 


Therefore,  I have  in  the  large  n limit  (s  ~ 1/n  •£  1) 


£|  (/*.»)  = 


C.(M  ^n^A  + v/5 
-Ct(p,*)A  t/s(l-A) 


(B.6) 


Here  ) = e ’ After  the 


Z,(M,Ar/2)  = Exp(SW)ErM 


(B.7) 


Combining  Eq.  B.l  and  Eq.  B.7,  I obtain  the  partition  function  Z(p,Af)  as 
shown  in  Eq.  2.9  in  chapter  2.  In  order  to  calculate  the  relative  weight  of  the 
configurations  with  k crossings,  I have  to  perform  the  inverse  Laplace  transform 
of  Z(ft).  Since  Z(/i)  behaves  like  1/p3,  the  relative  weight  of  configuration  with 
k crossings  is  ~ k/N. 


APPENDIX  C 

CALCULATION  OF  THE  CRITICAL  VELOCITY 


Here  I describe  the  evolution  of  two  small  hernias  until  the  disappearance  of 
one  of  them.  I calculate  the  average  lifetime  t of  a hernia  from  its  birth  until  its 
disappearance  and  average  it  over  all  possible  asymmetries  -1  < A < 1. 

I consider  the  equation  of  motion  for  two  hernias  with  length  L,  and  Lj(Fig. 
3.3).  I assume  stable  mass  input  into  each  segment  with  corresponding  rates  I\ 
and  I2  (in  dimensionless  variables  /,  = vl/va).  The  total  length  of  both  hernias  is 
Lt  = 2(Li  -f /.?)  and  the  change  of  Lt  is  proportional  to  the  total  input  /,  = Ix  +I2: 


dL,/dt  = /,v„ 


(C.l) 


If  the  influx  is  independent  of  time,  It  = constant  , then  Lt(t)  — /< Uo ( f + t0), 
where  £0  = Lt(t  = 0)/Itvo. 

The  equation  of  motion  for  the  asymmetry  factor  A = {Li  - Lt)/(Li  + is) 
can  be  obtained  from  Eqs.  3.2  and  3.3: 


where  the  influx  difference  is  defined  as  /.  = I\  — 1 2.  This  equation  can  be 
integrated  directly: 


(0.3) 


where  the  exponent  c = /,/( 2 - /,)  and  A,  is  the  initial  value  of  the  asymmetry. 
The  Eq.  C.3  can  be  solved  for  the  time  t it  takes  to  reach  an  asymmetry  factor 
A. 


When  one  of  the  hernias  disappears,  it  means  that  cither  L\  = 0 or  Li  = 0.  For 
definiteness,  I assume  that  the  second  hernia  disappears;  so  at  this  moment,  the 
asymmetry  A = 1.  Eq.  C.4  with  A = 1 gives  the  lifetime  for  this  hernia  as  a 
function  of  initial  asymmetry  A0.  I find  the  average  lifetime  by  averaging  t over 
all  values  of  initial  asymmetry  A0: 


(C.5) 


Note  that  the  integrand  is  divergent  in  the  vicinity  of  .-1„  - — /_/( 2 — /,).  For 
c < l(/(  < 1)  the  integration  over  this  singularity  is  convergent,  so  the  lifetime 
t is  finite.  It  means  that  the  newly  formed  hernias  disappear  after  some  finite 
time.  On  the  other  hand,  if  c > l(/f  > 1)  the  integral  is  divergent  which  means 
that  the  newly  formed  hernias  are  statistically  stable  (as  long  as  the  influx  stays 
constant  and  the  growing  hernias  do  not  attempt  to  split).  The  critical  velocity 
is  obtained  when  c = 1 with  symmetric  input,  /i  = h = 1/2-  This  results  in  the 
critical  velocity  vc  = v„/2. 
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